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O ■ Abstract 

o 

. We present a detailed examination of the variational principle for metric 

general relativity as applied to a "quasilocal" spacetime region M. (that is, 
a region that is both spatially and temporally bounded). Our analysis relies 
on the Hamiltonian formulation of general relativity, and thereby assumes a 
foliation of M into spacelike hypersurfaces S. We allow for near complete 
generality in the choice of foliation. Using a field-theoretic generalization of 
Hamilton-Jacobi theory, we define the quasilocal stress-energy-momentum of 
the gravitational field by varying the action with respect to the metric on the 
boundary dM. The gravitational stress-energy- momentum is defined for a 
two-surface B spanned by a spacelike hypersurface in spacetime. We exam- 
ine the behavior of the gravitational stress-energy-momentum under boosts of 
the spanning hypersurface. The boost relations are derived from the geomet- 
rical and invariance properties of the gravitational action and Hamiltonian. 
Finally, we present several new examples of quasilocal energy-momentum, in- 
cluding a novel discussion of quasilocal energy-momentum in the large-sphere 
limit towards spatial infinity. 
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I. INTRODUCTION 



Beginning with the earliest days of general relativity and continuing to the present, 
relativists have actively sought to define gravitational stress-energy-momentum (sem) from 
a variational principle. The motivation to do so is readily apparent, sem, and energy in 
particular, plays a central role in most branches of physics. In this paper we discuss a 
relatively new approach (see for example 1 Refs. [|IH25|]) to the problem, which we refer to 
here as the canonical quasilocal formalism (cqf). The CQF is based upon a field-theoretic 
generalization of Hamilton-Jacobi theory. We present many results new to the CQF and, in 
the process, recover the recent results from Refs. ]T2]JTE[ . 



Over the last thirty years, research has yielded a more-or-less satisfactory understanding 
of total energy-momentum for asymptotically flat spacetimes and asymptotically anti-de 
Sitter spacetimes. However, as no physical system is ever truly isolated, these asymptotic 
conditions — however useful — are ultimately unphysical, theoretical idealizations. In any 
case, practical numerical calculations are always restricted to a spatially finite region. For 
this reason and others (see the next paragraph), recent efforts have turned to the issue of 
defining sem quasilocally, that is to say, associating gravitational SEM with spatially bounded 
regions. (Indeed, one direct application of the formalism we present here is an approach to 
numerical outer boundary conditions for the gravitational field described in a forthcoming 
paper. pEfl ) 

As we will see, the CQF naturally leads to a definition of gravitational SEM that is 
quasilocal. We are motivated primarily by the desire to obtain physically meaningful and 
useful energy-like quantities that characterize the classical gravitational field in a bounded 
region. However, our original motivation for developing the CQF stemmed from a problem 
in semiclassical gravity, namely, understanding thermodynamical internal energy for black 
holes. The asymptotically-defined Arnowitt-Deser-Misner (adm) energy P?fl , for example, 
cannot serve as a useful internal energy because an infinite, self-gravitating system at finite 
temperature is thermodynamically unstable. Thus, the partition function can be defined 
only for systems with finite spatial extent, and this necessitates a quasilocal definition of 
energy (see Refs pB|-p(]|,PtPT|-|5^] and references therein). 

Before turning to the CQF, let us mention several approaches toward defining gravi- 
tational energy from a variational principle. The history of this problem is long, so an 
encompassing study would require a separate, extensive review. Here, we give only a brief 
summary of several of the historically important works. These works are based on a field- 



theoretic generalization of Noether's theorem |35 . 

Einstein was the first to derive gravitational SEM from an action principle. |36j By dis- 
carding a metric-dependent divergence term in the second-order covariant Hilbert action, 
he obtained a first-order action, the so-called IT action, that is the four-integral of a 
bulk Lagrangian quadratic in the Christoffel symbols. He then carried out a Noether- 
type analysis, and derived a canonical gravitational SEM pseudoteiasor and its corresponding 
super-potential. 2 Given what we've learned about the asymptotic structure of spacetime in 
the decades since this early work, it is remarkable how successful the Einstein definitions 
were. [37| Most of the key properties of spatial infinity (including decay of the metric and 
derivatives of the metric) are found in Einstein's original paper. The drawback of Einstein's 
approach is that the IT action is not fully diffeomorphism invariant (it is invariant modulo 
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boundary terms), and his gravitational SEM is coordinate dependent. At the quasilocal level 
there is no obvious general prescription for how one should choose coordinates. 

In the early 1960's M0ller discovered a new bulk action for general relativity that is 
similar to Einstein's IT action, but is quadratic in the tetrad connection (Ricci rotation) 
coefficients. |39| We might refer to it as the ujuj action. Like the Einstein Lagrangian, the 
M0ller bulk Lagrangian differs from the Hilbert Lagrangian by a pure divergence. (See 
also related work in Refs. [|0-H2].) Although the M0ller action is not fully invariant under 
"internal" transformations of the tetrad, it is diffeomorphism invariant and is therefore 
arguably preferable to the Einstein action as the starting point for a Noether-type analysis. 
Moreover, the resulting theory of SEM can be translated readily into the language of two- 
spinors, a powerful formalism which has led to numerous results in the quasilocal setting 
(see, for instance, the works of Szabados, Refs. [f43[-ff6f| and references therein). We point out, 
however, that in adopting the M0ller action, one is departing from pure metric relativity, 
Einstein's original theory. It is not clear at all that the SEM concepts derived in any such 
framework "pull back" to the metric phase space. 

By introducing a background metrical structure, one may isolate — in a coordinate in- 



dependent fashion — a purely metric divergence term in the Hilbert action. In Ref. |47 
Rosen discarded such a term, thereby obtaining another bulk action amenable to Noether 
techniques. An invocation of the Noether theorem in purely metric gravity, this approach 
towards defining gravitational SEM is close in spirit to Einstein's, and may be considered 
as a a refined version of his original analysis. However, the approach would seem limited 
in that there is not always a natural choice of background spacetime. Bicak, Lynden-Bell, 
Katz, and Petrov have developed and used an improved version of the approach in several 
recent papers (see Refs. f48[-50| and references therein) addressing, among other things, 



gravitational perturbations of cosmological solutions to the Einstein equations. 

We now turn to the canonical quasilocal formalism. Our analysis is based on the so- 
called "Trace-K" action pT],pLP], which differs from the standard Hilbert action by metric- 
dependent boundary terms. Its use leads to a purely metric formalism, as does the Einstein 
rr action. However, unlike the rr action, the Trace-K action 3 is manifestly invariant 
under coordinate transformations. Moreover, the Trace-K action does not depend on any 
background structures. 4 Since the Trace-K action does not stem from a bulk Lagrangian, it 
is not immediately clear how to apply the Noether theorem. But we can bypass a Noether 
analysis altogether, using instead the CQF which is based on Hamilton- Jacobi theory. We 
point out that our approach is intimately related to a body of work done by Kijowski and 
co-workers (see Ref. |52| and references therein). Kijowski's approach starts with novel 
and important ideas from symplectic theory ||53|| , and examines the relevant symplectic 
geometry in great detail. Our approach, on the other hand, starts with standard Hamilton- 
Jacobi theory and focuses primarily on the physical spacetime geometry. We stress that 
both approaches are merely different faces of a Hamiltonian analysis, and thus somewhat 
different from more traditional approaches based on Noether techniques. 

Consider a spatially and temporally bounded spacetime region M. with metric g^ u and 
boundary dAi. The boundary dM. of the region consists of a timelike element T (the 
meaning of the "bar" is explained below) and spacelike elements E' and Such a spacetime 
region is depicted in Fig. 1, but note that T need not be connected. We assume that the 
spacetime M. is foliated into spacelike hypersurfaces E, defined by t = const. Further, 



3 



FIG. 1. 



we require the boundary of each E leaf to lie in T. 5 The intersections of the leaves of the 
spacetime foliation X with the timelike boundary element T define a foliation of T into 
two-dimensional spacelike hypersurfaces B (which need not be connected). The region M. 
may itself be contained in some ambient spacetime. We note that the boundary B and its 
history T are simply submanifolds of spacetime and need not be physical barriers. 

The future-pointing unit normal to the t = const hypersurfaces is denoted w M , and the 
outward-pointing unit normal of T is denoted n M . The induced metric on S is hij and the 
induced metric on T is 7^. Because in general n^u^ ^ on T, the Eulerian observers of 
the B foliation of T, comoving with T, need not be at rest with respect to the S slices. 
This necessitates the "barred" and "unbarred" notation which keeps track of the two sets of 
Eulerian observers at the boundary: those comoving with T and at rest with respect to the 
B foliation of T (the "barred observers") and those at rest with respect to the £ foliation 
of M. (the "unbarred observers"). 

The four-velocities of the barred observers will be denoted u^ 1 . Similarly, at each point 
of T, we define n M as the unit outward pointing vector for the unbarred observers. That is, 
n M lies in S and is orthogonal to B. 6 Note that, by construction, u^n^ = and w M n M = 0. 
These vectors are related by the boost relations 

% = 7?V + l vu ^ > (1.1a) 
= 7?V + l vn v ■ P) 

where v is the boost velocity between the two sets of observers and 7 = (1 — v 2 )^ 1 ^ 2 . In 
Appendix A we present the details of the kinematical relationships needed for this paper. 

As mentioned, our analysis is carried out in the purely metric formulation of gravity and 
is based on the Trace K action, [|^ ,P,|^|T7| 

S[g] = — f d 4 x^^+ - F d 3 xVhK - - f d 3 xy/^Q - - [ B d 2 x^9 . (1.2) 
2k Jm k js' k Jr k Jb 1 
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Here, k is 8ir times Newton's constant. For simplicity we have omitted matter and cosmo- 
logical constant contributions to the action. The symbol J s , is shorthand for J s „ — J s ,, and 
K is the trace of the extrinsic curvature K^ u = —h^SJ a u v of the boundary elements E' and 
E". Similarly the function 6 is the trace of the extrinsic curvature = — a n v of the 
boundary element T . The action ( |1 . 2| ) includes contributions (first considered in Refs. [§,[§) 
from the "corners" B' = E' PI T and B" — E' D T, where <r is the determinant of the metric 
a a b on the corners. The velocity parameter 9 is defined by sinh# = —u^n^ = jv. 

The Trace-K action has the crucial property that its associated variational principle 
features fixation of the induced three metric 3 gij on dM.. 7 In particular, the lapse of proper 
time for an observer, comoving with T and at rest with the B foliation, is fixed as boundary 
data since this information is encoded in the fixed T three-metric. The value of the quasilocal 
energy surface density (at a given point on the observer's worldline) is defined through a HJ 
variation as minus the rate of change of the classical action with respect to an infinitesimal 
stretch (enacted at the given point) in the proper time separation between E' and E". Of 
course, the T three-metric specifies more than just the lapse of proper time between the 
initial and final slices — it contains information about all possible spacetime intervals on T. 
One is free to consider the changes in the classical action corresponding to arbitrary HJ 
variations in the T metric. The original Ref. [|I|] has demonstrated how this freedom leads 
not only to the energy surface density but also to surface densities for tangential momentum 
and spatial stress (both are pointwise-defined B tensors). 

In this paper, we extend the CQF analysis of Ref. |l| by considering changes in the classical 
action corresponding to HJ variations in E" (or E') boundary data. This leads to quasilocal 
surface densities for normal momentum, tangential momentum (which is equivalent to the 
previous definition), and temporal stress (this was also shown in Ref. [|12|)- Therefore, the 
quasilocal stress-energy-momentum consists of energy, normal momentum, and tangential 
momentum surface densities and spatial and temporal stress tensors. 

We also extend the CQF by considering "boost relations" between the quasilocal surface 
densities as defined by barred and unbarred observers. These SEM boost relations can be 
viewed as canonical realizations of the relations ( |1 . 1|) satisfied by the barred and unbarred 
observers' unit vectors. The observer dependence of the quasilocal SEM is best described from 
the following perspective. The various SEM quantities are defined as tensors on the spatial 
boundary B spanned by a spacelike hypersurface E. The boost relations characterize the 
behavior of these tensors under a boost of the spanning slice E; that is, they characterize 
the dependence of the quasilocal SEM on the choice of observers passing through B. As 
purely geometrical relations, the boost relations among energy, tangential momentum, and 
normal momentum surface densities have been noted elsewhere in the literature (see, for 
instance, Refs. fli3| , |5^1 ). Moreover, their particular role in the CQF has been pointed out in 
Refs. H5,12| . Here we present a unique derivation of these relations, demonstrating that their 
geometrical content is already encoded in the gravitational Hamiltonian. 

Finally, we present several new examples of quasilocal energy-momentum, including an 
analysis of cylindrical gravitational waves. We also include a novel discussion of quasilocal 
energy-momentum in the large-sphere limit towards spatial infinity. Agreement between the 
Trautman-Bondi-Sachs total graviational energy-momentum and the notion of quasilocal 



energy-momentum arising in the CQF has been considered elsewhere. [16 



In Sec. II we prepare for the Hamilton-Jacobi variation of the Trace-K action by con- 
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sidering the general variation of the Hilbert action. Of particular importance are the corner 
terms that arise at the intersections of T with £' and These terms have appeared 
previously in the literature. P,|3| ]T2| . |52|| However, to our knowledge, Sec. II contains the first 
completely geometrical derivation of the result (although the same result was obtained ex- 
plicitly via another method in Ref. ||52|| ). In Sec. Ill we apply the CQF to the Trace-K action 
and derive the quasilocal sem. In the process, we obtain the boost relations among the en- 
ergy and momentum surface densities and spatial and temporal stress tensors as defined by 
the boosted and unboosted observers. We also discuss the notion of boost invariants which 
allow for the construction of several mass definitions which have appeared in the relativity 
literature. Section IV contains a derivation of the Hamiltonian form of the action and its 
variation. We then derive the boost relations for the energy and momentum surface densities 
by boosting the gravitational Hamiltonian. Finally, Sec. V and Sec. VI deal with concrete 
examples and applications, and in these sections we chose units such that the constant k 
appearing in Eq. ( |1.2| ) is simply 8tt. In Sec. V we first deal with the issue of zero-points 
for the quasilocal densities, and the relationship between these zero-points and a freedom 
always present in any variational principle. The freedom concerns appending to the action, 
here the gravitational Trace-K action, an arbitrary functional of the fixed boundary data. 
Having (at least partially) dealt with this issue, we then write down quasilocal energy and 
momentum expressions for a variety of exact solutions to the Einstein equations; and in 
Sec. VI we apply our formalism to spacetimes with are asymptotically flat in spacelike di- 
rections, in the process making some novel observations about total gravitational energy at 
spatial infinity. 

Appendix A contains several key kinematical results that are used throughout the paper. 
Appendix B is devoted to the derivation of certain curvature splittings needed for the analysis 
in Sec. III. Finally, in Appendix C, we show that the rate of change of the boost parameter 
equals the normal gradient of the lapse function defining the boost. This is needed for the 
analysis in Sec. IV. 



II. VARIATION OF THE HILBERT ACTION 



In this section we consider the standard Hilbert action, [54 



S H [g\ 



2k Jm 



(2.1) 



and its associated variational principle as applied to a bounded spacetime region A4, a 
careful analysis of which is crucial for the entire discussion. Such an analysis is, of course, 
not new |51],^|fjrr^ri^j5^] ; however, as we do give a new version of a nontrivial calculation 
of fundamental importance, we believe the details belong up front and not relegated to an 
appendix. 

The relevant geometry of the various foliations of A4 is described in the Introduction 
and Appendix A. We examine the variation 8S H of the action induced by an infinitesimal 
variation dg^ in the metric tensor and derive the following result: pT],p|,|3|, p^?| , p!8| , |52 
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5S H = -— [ d 4 xJ=g~G^5g, lu - F d^xh^P^ - [ d 3 xf ii 8U ij + - F d 2 x^56 . 

2k J M JT,' JT K JB' 

(2.2) 

In this expression G^ v is the Einstein tensor, 9 is the velocity parameter described earlier, 
and 

pv = ^ (Rh l] - K ij ) (2.3a) 



2k 

Y[ ij = -^y^- ~ (lib) 



are respectively the £ and T gravitational momenta. 



In writing the Hilbert action ( |2.1| ) and its variation ( |2.2| ), we do not necessarily assume 
that the spacetime boundary elements X", and T have smooth embeddings in M.. That 
is, the unit normal n M of T and the unit normals of £' and S" need not be continuous 
vector fields. For example, the timelike boundary element T can contain a "kink" , a spacelike 
two-surface at which changes discontinuously. In that case the T boundary term in 5S H 
contains a contribution from the kink. The form of the contribution is discussed in Ref . || . 



A. Preliminary results and a lemma 



To begin, let us collect a few results concerning the variation ST \ u of the affine connec- 
tion induced by an infinitesimal 5g fiu . With these we prove a lemma of particular use when 
examining the variation 8S H of the action. 8 First, expansion of the identity 5(V\g^ u ) = 
leads directly to the result 

2<ffV )A = V x 5g^ . (2.4) 



Second, the well-known formula I531 cJJLj, = V\5T X au — V^r^A for the induced variation of 



the Ricci tensor implies that the contracted variation g tJ " u 5$t, lll/ is a pure spacetime divergence. 
Indeed, writing g^d^R.^ = V^ M , we find that 

V fM = 25T^ u] u . (2.5) 

Finally, consider a metric-dependent covector (i. e. du^ need not vanish) and the space- 
time divergence V^uj 11 constructed from it. The variation of this divergence is 

= -i^y + - \{V»u»)5g, v - \V»{u v 8 9lxv ) . (2.6) 

To obtain this result, expand the variation 5(g (ll 'V \tUJ u ), insert the identity tfV^cjj, = V^5lj u — 
u\5T x llu , and then use Eqs. ( [2.4|) and (|2.5| ). 

We now prove the following. 
LEMMA : Consider a unit hypersurf ace-orthogonal vector field, say with normalization 
u^u^ = e (with e = ±1). 9 Also consider the induced metric h^ u = g^ u — eu^u u on the 
hypersurfaces to which w M is orthogonal, as well as the extrinsic curvature tensor K^ IV = 
—h^VaUf- The variation 5K of the mean curvature K = satisfies the equation 
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25K = - K^Sh^ - D a (h«6u») , (2.7) 

where D a is the covariant derivative operator compatible with h^. 

Our proof of the lemma makes use of the following two identities: (i) W = —K^ v + 
eu^u x V ' \u v and (ii) 5u fM = ^eu IM u a u l3 6g al 3. Identity (i) follows directly from the definition 
of K^ v and the spacetime expression for the induced metric h^ u . We verify (ii) by writing 

= eNV^t, where the coordinate t labels the hypersurfaces to which w M is orthogonal and 
N = (eVfj,tg^ u Vyt) _1//2 is the lapse function. As the first step towards proving the lemma, 
we rewrite Eq. (|2.6|) with in place of u^, make substitutions with the identities (i) and 
(ii), and do a bit of algebra in order to obtain 

25 K = UpV - K^Sg^ - e(V x u x )u>"u v 5g^ + W X V \{u v 5g^) . (2.8) 

Now, by (i) V M u M = h^^^Uu', and, therefore, we may collect the last two terms on the 
right-hand side in Eq. ( |2.8|) , thereby arriving at 

25K = n^V - K^8 9fll/ + /i Aft V \{hfy8g^) . (2.9) 

Finally, since K^ u is purely spatial, K^Sg^ = K^Sh^. Moreover, with identity (ii) we 
can show that h^u v 5g^ u = —h^du^. Substitution of these results along with the definition 
of D a into Eq. ( |2.9| ) completes the proof. 

As we have been careful to allow for the case e = 1, our proof of the lemma establishes 

250 = - e""#v - V a {%8^) (2.10) 

as a corollary. Here, T> a is the covariant derivative compatible with the T metric 7^. 



B. Variation of the action 



Our goal now is to obtain the expression (|2.2|) for the variation 5S H of the Hilbert action. 
Simple manipulations show that the variation of the action is 

5S H = ±- [ d 4 x^g- (-G^Sg^ + V„V) , (2.11) 
2k Jm 

where has been defined in Eq. ( |2.5|) . Focus attention on the divergence term in Eq. (|2.11 ), 
namely, 



Via Stokes' theorem, 



(5S H ) dM = / d A x^-gV^. (2.12) 
2k Jm 

JM JdM 



we can express {SS h )qm as a pure boundary term. In Eq. ( |2.13| ) we have used abstract 
notation for both the spacetime volume form (e) IJV x K = e^vx*. and the three-form (V -e) uXK = 
y^^iivXn- The integral on the right-hand side of Eq. ( j2.13| ) can be written as a sum of separate 
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integrals over each element of the boundary dM = E' U E" U T . Indeed, no corner-term 
contributions (i. e. two-surface integrals over B' and B") can arise at this stage, because 
both B' and B" are sets of measure zero with respect to dM and the integrand (V ■ e) 
is continuous (as both V 11 and e^ u \ K are continuous). Now, standard convention fixes the 
orientation of the boundary dM by choosing the outward-pointing normal n M to dM as 
embedded in M] that is to say, the alternating tensor on dM is taken to be e v \ K = n^e^xn- 
Subject to this convention, one finds in general that 

/ (V-e)=[ dPxJfgWn.V^ . (2.14) 

JdM JdM V 

In this expression \ 3 g\ is the (absolute value of) the determinant of the induced metric on 
dM and e = n M n M is a sign factor which is either 1 or —1 depending on the boundary element. 
Notice that en M is the covector dual to the outward-pointing normal n M . Therefore, for the 
case at hand with dM = E'IJE"IJT, we expand the right-hand side of Eq. Q2.14|) , and 
obtain 

(SS H ) aM = ~ C d 3 xVhu^ + /_ d 3 x^/n^ (2.15) 



as the promised boundary expression for the divergence term (1212]) . 



Next, combining the lemma ( |2.7| ) and its corollary ( |2.10 ) with ( 2.15 ), we write the 
divergence term as follows: 

1 r B " 

(SS H ) dM = -— J d 2 xy/a (n fl 5u tl + u^) 

-— F d 3 xVh{25K + K ij 5h ij ) + — f d 3 x^/^f(25Q + 0^5%) . (2.16) 

In this expression the total derivative terms in the boundary element integrals have been 
expressed as integrals over the corners B' and B" (again via Stokes' theorem, but now in one 
lower dimension). Further, for the E', E" and T integrals in Eq. ( 2.16|) we have chosen to 



use indices adapted to the boundary elements rather than general spacetime indices. Now, 
with the momenta (|2.3|) used in Eq. (|2.16|), the term (SS h )qm becomes 



1 r B " r 

(SS H ) dM = / <Px^(nJu* + uJn>>)- / d 3 xh ij SP lJ - / d 3 xj ii SU v . (2.17) 

2k Jb> Jt,' Jt 

Our remaining task is to simplify the integrals over B' and B" . To achieve this, recall the 
boost relations 

v? = iii" + iryn" (2.18a) 

7n M + vju^ (Hi b) 



and their inverses, derived in terms of a double foliation of spacetime in Appendix A. These 
can be used to write the integrand of the corner integrals as 

n^Su 1 * + u^n^ 1 = (n M /7 - vu li )6(W/'y - vn^) + (u M /7 + vn tl )5(n^/'y + vu^) . (2.19) 
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On the right-hand side of this equation, the terms proportional to u^dn^ vanish. This 
follows from the identity u^dn^ = —n^Su^ and the fact that is hypersurface-orthogonal. 
[Thus, as seen in identity (ii) after the Eq. (|2.7|) 5u^ is proportional to u^] Likewise, we 
have n^Su^ = 0, since n M is hypersurface-orthogonal. After a bit of straightforward algebra, 
Eq. ( |2.19| ) simplifies to n^du^ + u^n^ = —2^ 2 bv. Therefore, we may now rewrite Eq. ( [2.17|) 
as 

1 r B " /-s" r 

(SS H ) dM = - / d 2 x^-f 2 5v - / d 3 xhii5P v - / d 3 x^ ij SU v . (2.20) 
k jb 1 Jy,' Jr 

Combination of this result with Eq. ( |2.11| ) and the definition tanh# = v of the boost pa- 
rameter yields the desired expression ( |2.2j ). 



C. Boundary terms and the diffeomorphism invariance of the Hilbert action 

The Hilbert action ( |2.1| ) is diffeomorphism invariant. That is, the action is unchanged 
if the variations in the fields are given by the Lie derivative along a vector field £ M that is 
tangent to the boundaries: 



H = ±-( d A x £s(V=9®) = 7T / d±xV„(eV=9®) 
2k Jm n 2k Jm 

[ d 3 xVhuJ^+ — [ (Pxyf^n^ 11 ®. 
Jy,' 2k Jr 



2k 

= 0. (2.21) 

Here, we use = on T, and = on E' and E". These imply = = on 
B' and B" . 

Since 5Sh = when 5 is given by the Lie derivative, our main result Eq. ( p.2[ ) implies 

= 6S H = ~- I d A x^G^Vaiu + - l B d 2 x^ 1 2 Cd a v 
k Jm k Jb> 

-J d 3 xh ij £ ( P v - J_d 3 xf ij £ i U v . (2.22) 

Now use the identity 

h, r (:J" J = I J> - l" J tJi,, = Di(PC) - 2P ij D i £ j 

= A(Pe-2P%)+2(AP i % , ' (2.23) 

where is the pullback of £ M to £' or E". Note that the _Dj(P£*) term will vanish when 
integrated to the corners. Thus, the £' and E" terms in ( |2.22| ) become 

rY" rY" rB" 

- J d 3 x h i:j £^P ij = -2 J f d 3 x (/;, /"')(.; + 2 J d 2 x^/a ni P ij ^/Vh . (2.24) 
Similarly, we find 

- J d 3 x "iytAV = -2 J d 3 x (ViW 3 )^ -lf Bi d 2 x^u i Wijl4^( , (2.25) 
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whence Eq. (|2.22 ) now becomes 



= 5S H = -- I <PxV=jj{V + f B d 2 x^ (U a ~ JaC + l 2 Cd a v/K) 

k Jm j b> v ' 

+ d 3 x(Vhu fl G^^/K-2(D i P i %) 

- j_d 3 x(^n^ v £ v /K + 2{V i Il i %) , (2.26) 



where we have used integration by parts on the volume (M) integral term. Also, we have 
used the definitions j a = — 2o a i IP 3 ' Uj / -y/ — 7 and j a = —2a ai P^nj/\/h. 

We now use the well-known result that the gravitational field contributions to the bound- 
ary momentum constraints satisfy 



-2D 3 P( 



Tti 



n, = -2P,ni 



-VhupG^hvi/K 



[2.27) 
[2.28) 



Therefore the last two integrals in Eq. ( |2.26|) vanish. Since the result ( 2.26 ) must hold for 
all ^ that are tangent to the boundary, we conclude that 



V,G^ = , 

Ja=Ja 9 a 9 

K 



(2.29a) 
fljgg b) 



where 6 is the velocity parameter, v = tanh(6 l ). Equation (|2.29| a) is, of course, the contracted 
Bianchi identity. Equation ( 2.29| b) is an identity as well. In fact, as we will see in the 
next section, j a and j a are the tangential momentum densities for the barred and unbarred 
observers, and the identity ( |2.29| b) expresses the boost relationship between these quantities. 

Note that this analysis can be applied to the Trace-K action as well. Indeed, any action 
that is diffeomorphism invariant and differs from the Hilbert action by boundary terms can 
be used. The reason is that the Lie variation of a boundary term will always integrate to 
the corners, and then vanish since £ M is tangent to the corner. 



III. QUASILOCAL STRESS-ENERGY-MOMENTUM AND BOOST RELATIONS 

A. Quasilocal quantities 



Using our main result (|2.2) for the variation of the Hilbert action, one can easily show 



that the variation of the Trace-K action (1.2) has the following boundary terms: 



(5S) dM = r d 3 xP ij 8h i:j + I d 3 xU^8j i:j ( B d 2 x65^ . (3.1) 

JE' JT K JB' 

Notice that the Trace-K action features solely fixation of the induced metric on the boundary 
dJ\4. We now wish to express the T boundary term in SS in terms of the geometry of the 
£ slices. Start with the {SS)f contribution to the variation, that is 
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T 



d 3 x U ij S^ 



(3.2) 



With an ADM splitting of the side boundary metric 7^ into a lapse function N, a shift 
vector V a , and a spatial metric a ao (see Appendix A), we find 



Sjij = —=-UiUj8N - T -a a{i u j) 5V a + a^aida, 



ab 



N J N 
With this splitting of the T metric we then obtain 



d 3 xU ij 6f ij = -- J_d 3 x^\ (9 + u0 j Uj)5N + (aaiQ^Uj^V 



(3.3) 



r 



.06 



a 



(7 



<?ab 



(3.4) 



Now, in order to achieve our goal of expressing (5S)f in terms of £ geometry, we must find 
a "splitting" of the T extrinsic curvature tensor Ojj. The desired expression 







ik^ + jvi^ + (n ■ a)u^u u + 2a*u v} (K a pn p 



is derived in Appendix B. In this expression, we have used the definitions 



(3.5) 
(3.6a) 



The unit normals w M , n M associated with the hypersurfaces £ are related to the unit normals 
u^, n M associated with T as in Eqs. fl2.18|) . Again, our conventions are that barred observers 
are comoving with the boundary T while the unbarred ones are at rest in the £ hypersurfaces. 
Also in Eq. (3J5), a M = u u 'V l ,u >1 denotes the acceleration of the barred observers, and K a p 
denotes the extrinsic curvature of the £ slices. Putting these results together, we have 



d 3 x \\ ij (S', / = f d 3 x^a~\ [yk + jv£)6N - KK ljn j - d a 6}5V a 



N 



~i{k ab - ka ab ) + n-da ab + 1V {t b - £a ab ) 



— — ^ab 



ab 



r ab\ 



5a, 



ab 



(3.7) 



for the T term in the variation of the action. 

The contribution to 5S from the top and bottom caps (£" and £') is 

/ d 3 xP lJ 5hij = — d 3 xVh(Kh ij - K ij )6h i:i . 



(3.8) 



The induced metric can be split into a "radial" lapse function M, shift vector W a , and 
slice metric a a b (see Appendix A). The variation in hij is then 



Sh i:j = — + ■^a a(l n j) 5W a + a^Sa, 



M 



M 



(3.9) 



from which we obtain 
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f 

Jt< 



d 3 xP lJ 8h, 



K JT 



T" 



d 3 xy/a\ (K - mK^n^SM - (a ai K ij nj )8W a 



+^-{Ka ab - atK ij a))5a ab \ (3.10) 



Now use the splitting 



K, 



+ {u ■ b)n^n v + 2a? n v) K a pnP 



(3.11) 



from Appendix B, where = n u 'V u n fl . These results yield 



Jt< 



d 3 xP ld 8h 1 



K JT 



T" 



d 3 x^a~\ i8M - a a K ijn j 8W a 



M 



{r -io ab )-u-bo 



ab 



(3.12) 



for the top and bottom-cap terms in the variation of the action. 

The result (|3.7| ) allows us to define the quasilocal densities associated with the two- 
surfaces B as seen by the "barred" observers: 



K£ 



KS ab = 



k 8S\ 



T 



fa 8N 
k 5S\f 

8V a ~ 
2k 8S\ f 
N^fa 8a ab 



7& + jvi 



a\K i5 v? - d a 9 , 



7(fc a " -ka ab )+n-aa ab +jv(i 



ab 



f,ab 



ab 



(3.13a) 
(|3TTB|b) 
( gT3| c) 



These are the quasilocal energy density, tangential momentum density, and spatial stress, 
respectively. The notation 8S\f refers to a HJ variation of the Trace-K action S, with respect 
to the T metric components N, V a , and a a b- These definitions hold for each leaf of the T 
foliation B, but our attention will be focused primarily on the corner B" . Likewise, the 
result ( |3.12|) allows us to define quasilocal densities as seen by the "unbarred" observers: 



k 8S\ : 



a 5M 
k 8S\ E „ 



Kt 



ab 



fa 8W a 
2k 5S , | s „ 
M^fa 8a ah 



-a^Ktj = -I 



a a K ijn i 



Ka ab - a?K ij a b = -[t b - ia ab ) +u-ba 



ab 



(3.14a) 
flgTg b) 



These are the quasilocal normal momentum density, tangential momentum density, and 
temporal stress, respectively. The notation 8S\^„ refers to a HJ variation of the Trace-K 
action S with respect to the E" metric components M, W a , and a ao . These definitions hold 
for each slice of the "radial" foliation of E", but again we focus attention on the corner B" . 

Clearly the definitions (|3.13| ) and (|3.14|) are applicable to any closed two-dimensional 
surface B embedded in a spacetime that satisfies the Einstein equations — we simply arrange 
to have the top corner B" of the manifold M. coincide with the given surface B and apply 
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the definitions. The surface B can be pierced by various fleets of observers, for example, 
barred and unbarred observers. Different observers who are boosted relative to one another 
will see different quasilocal densities for the same surface B. With this in mind and to put 
the set ( |3.14|) on an equal footing with the set (|3.13|) , we define additional barred densities 



nj h = -a ij Kij = -ji - yvk , (3.13d) 
Kj a = aiKijfv 1 = a l a K tJ n> - d a 9 , ( [Ogl e) 
nt ab ee Ka ab - a^K ij a b = - 7 (r b - £a ab ) +u-ba ab - ^v{k ab - ka ab ) . fl£I3| f) 

Note that these expressions are defined in terms of a slice £ [with intrinsic and extrinsic 
geometry (hij, K^)} which meets the T boundary orthogonally. Observers comoving with T 
are at rest with respect to £. 10 It is not difficult to see that in terms of the (hy, K l i) geometry, 
Jh, J a, and i ab have exactly the same forms as do j a , and t ab in terms of (hij, K^) geometry. 
In the rightmost expressions we have expressed Jj_, j a , and i ab in terms of £ geometry, using 



a "splitting" similar to the one given in Eq. (|3.5|) but this time expressing the spacetime 
representation K^ u of the £ extrinsic curvature tensor in terms of £ geometry. The relevant 
splitting is found in Appendix B. Finally, note that expressions ( |3.13[ b) and ( |3.13| e) agree. 



In Eqs. ( |3.13| ), the quasilocal densities for the barred observers are expressed in terms of 
the geometry and foliation defined by the unbarred observers. Alternatively, those densities 
for the barred observers can be expressed in terms of the geometry and foliation defined by 
the barred observers themselves. This is achieved by keeping the boundary dM. fixed in 
a neighborhood of B", and tilting the £ slices until the unbarred observers coincide with 
the barred observers. In other words £ slices become £ slices. The boost velocity v then 
vanishes, and Eqs. (|3.13Q become 

K£ — k , (3.15a) 
Kj h = -I , (l315j b) 

KJ a = (JaKijU 3 , ( |3.15| c) 



K s ab = k ab - ka ab + (n ■ a)a ab , (|3TT5| d) 
Kt ab = -l ab + to ab + {u ■ b)a ab . ( gT5] e) 



Here, the barred quantities k ab , Kij, etc. refer to the surface B" embedded in the top cap 



S". The results (|3 . 1 5 ) extend the definitions given in the original QLE paper [HJ 11 to include 
the normal momentum density and the temporal stress tensor i ab . Of course, we can view 
the limit v — > of Eqs. ( |3.13|) in another way: consider the unbarred observers (S slices) 
as unchanged, and the boundary dA4 at the corner B" as "unboosted" until the barred 
observers coincide with the unbarred observers. Then we obtain the relationships (ft.lSQ , 



but without the bars. That is, we find that the energy surface density for the unbarred 
observers is ke = k, with similar expressions for the momentum densities and stress tensors. 

Before continuing with the main line of reasoning, let us discuss the physical significance 
of the normal and tangential momentum densities. The normal momentum density can be 
written as njy- = er^V pU v = —a^K^ = n l n^{Kij — Khij), and the tangential momentum 
density can be written as nj a = a % a n 3 {Kij — Khij). These quantities are the normal and 
tangential components of the (total) momentum surface density nf = rij{K 1 ^ — Kh^), 
which can be written in terms of the gravitational momentum as f = —2P l ^nj/\/h. We 
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now remark that the analysis presented in this paper can be easily generalized to include 
matter fields. For the case of nonderivative coupling (in which the matter action does not 
contain derivatives of the metric) the basic definitions ( |3.15| ) are unchanged. By including 
matter fields in the definition of the system we find that f = j^n 1 + ] a <2 a% is related to the 
matter momentum in the following way. Consider the momentum constraint 

= -2DiP ij + VhunTM (3 16 ) 

for the hypersurfaces E, where —u^T^ is the proper matter momentum density in the jth 
direction. Assume that there exists a Killing vector field on space E. It is straightforward 
to show that the total matter momentum along £ l is 

- f SxVhu^^i = -2 f d 3 xDiP ij ^ = f d 2 Xv ^3iC , (3.17) 
where B = <9E. This shows that f represents a surface density for the matter momentum. 



B. Boost relations 



We now return to Eqs. ( 3.13 ), expressing the quasilocal densities for the barred ob 
servers. In terms of the unboosted quasilocal densities [the unbarred versions of Eqs. ( |3.15| )1 
Eqs. ( |3.13| a,d) become 



e = 1 e- >yvj h , 



(3.18a) 
iHb) 



relating the quasilocal energy density and the normal momentum density for barred and 
unbarred observers. We also obtain 



Ja= 3a~ d a 0/ « , 

for the tangential momentum density. Finally, we have 



-ab 



— [n ■ a)a ab '/ k 



- ab 



- (n ■ a)a ab /K 



7 

i ab - (u ■ b)a ab /n = 7 t ab - (u ■ b)a ab /K 



— jv 

— -yv 



i_ab 



-,ab 



- (u ■ b)a ab /K 

- (n ■ a)o ab /n 



(3.19) 
(3.20a) 



for the boost relation between the spatial and temporal stress tensors. This later relation 
can be rewritten using the results 



n ■ a = 7n • a — ^vu ■ b + u ■ V8 
u -b = 7« • b — jvn ■ a — n ■ V6 



from Appendix B. We thus obtain 



-ab 
Tab 



.,ab 



7S- - ^vt ab + (u ■ V6/k)ct 
lt ab - 1V s ab - (n ■ V6K)a ab 



ab 



(3.21a) 



(3.22a) 
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for the boost relation satisfied by s ab and t ab . Finally, let us define the spatial shear r] ab = 
s ab — sa ab /2 as the trace-free part of the spatial stress s ab , and the temporal shear ( ab = 
t ab — ta ab /2 as the trace-free part of the temporal stress t ab . From the unbarred version of 
Eq. (13.151) , we have 

K7] ab = k ab _ ka ab/ 2 ^ (3.23a) 

K( ab = ~t b + £a ab /2 . flpg b) 



The results fl3.20|) , or equivalently (|3.22|) , yield 



rf> = ir] ab - ~fv( ab , (3.24a) 
( ab = 7 C ab - ivr] ab , flpg b) 



for the boost relation between rf b and ( ab 



C. Boost Invariants 



The results (|3.18| a,b) show that the energy surface density and normal momentum density 
behave under local boosts like the time and space components of an energy-momentum 



vector, namely, eu^+j^n^ 
defined by 



eu^+ji-n^. Clearly, the squared length of the vector eu^+ji-n^, 



M 2 /k 2 



(3.25) 



is invariant under boosts. We do not claim that M 2 is in all cases positive. However, if it is, 
then M/k (defined via the negative square root ]13[) is equal to e for a fleet of observers 
who pass through B in such a way that f\- = —£/k = 0; that is, such that B is a maximal 
slice of T (if such a slice exists). This defines locally, at each point of B, a rest frame for 
the system. Moreover, the parameter associated with the local boost between an arbitrary 
frame and the rest frame can be computed via 



m 



= §log 



(3.26) 



Indeed, using the relations inverse to those given in Eqs. (|3.18| a,b) along with the rest-frame 
condition jy. = 0, we find that 



1 +v 



which is immediately recognized as the logarithmic representation of tanh 



(3.27) 



Note that 



Eq. ( |3.26| ) demonstrates that the two-surface data encodes the rest frame direc- 

tion. This fact features prominently in Kuchaf's examination of the geometrodynamics of 
Schwarzschild black-holes. [|5l||T3|j 

Equation (|3.19|) expresses the change in the tangential momentum surface density j a 
under a boost. Evidently the curl of j a , 
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Fab = d a J b - d b J a , 



(3.28) 



is invariant under boosts. Also note that j a itself is invariant under boosts that are constant 
on B. 

We now turn to the spatial and temporal shear. The boost relations ( |3.24| a,b) show that 
the shear tensors transform like the components of a (two-dimensional, traceless, symmetric 
matrix valued) spacetime vector H ab ^ = k^tj^u^+C The shear tensors can be combined 
to form boost invariants, such as 

(Mi) 2 = H ab »H ab , = K 2 ( V ab Vab - CXab) , (3.29a) 
(M 2 ) 4 = H^HfH^H^ - (H ab »H ab ,) 2 = 2K* V ab ( Vab ( cd - U?lcd)C, cd , O) 
(M 3 ) 2 = e ab ^H c a ,H bcu = 2K 2 e ab V c a ( bc , (g) 

where e ab = u fl n l 'e^ uab is the alternating tensor on B. Again, we do not claim that (Mi) 2 , 
(M 2 ) 4 , and (M 3 ) 2 are positive. 

With the invariants (Mi) 2 , (M2) 4 , and (M3) 2 we can build several different mass defini- 
tions that have appeared in the literature. |12j In terms of the expansions p and p of the 
null normals (u^ ± n^) / \[2 to B (spin coefficients in the Newmann- Penrose formalism [|57][ ), 
we have that M 2 = 8pp. Hence the familiar Hawking mass [0] may be expressed as 



M Ha wking = ~Jtzt- I <l'--r 



K V 167T JB 



[ d 2 x^\-{M 2 /2) +TZ] , (3.30) 

J B 



where 1Z is the Ricci scalar and A the area of B. The prefactor in front of the integral 
ensures that the overall expression has units of inverse length (i. e. energy in geometrical 
units). 

The following is a geometric identity relating the M. Riemann tensor 3t a i3fiu with the 
two-surface data of B: |43| 



a^a x ^ XuK = k ab k ab - k 2 - i ab t h + £ 2 + n= {M x y - (M 2 /2) + K , (3.31) 
where here the B two-metric a^ v = g^ v — n^n v + u^u v serves as a projection operator. 



Hay ward's quasilocal mass [59 



is 



M Hayward = - J — I d 2 x (vWV Ak 3^ A(7K + a^u, nf[u, n] v ) , (3.32) 

where [u, n]^ is the vector-field commutator between the B normals. One may verify that 
the last term in Hayward's mass is boost invariant, although it would not seem expressible 
solely in terms of the two-surface data of B. Striking this term from the integrand one 
obtains an energy expression which has proved useful in asymptotic investigations. [70 

D. Second Fundamental Form of B in M. 

The second fundamental form (extrinsic curvature) for a spacelike two-dimensional sur- 
face B embedded in four-dimensional spacetime M. is defined by 
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#a/ = <r>* V 7 a£ . (3.33) 

Here, as always, a a p is the induced metric on B and V is the covariant derivative in Ai. 
With the representation o a p = g a p + u a up — n a rip, the second fundamental form becomes 

Ha^ = ala^V^us - n^V 7 n 5 ) . (3.34) 

From this result it follows that u^H^ = £ a p is the extrinsic curvature of B as a surface 
embedded in T, where T is the three-dimensional spacetime orthogonal to n M . It also follows 
that n^Ha^ = k a p is the extrinsic curvature of B as a surface embedded in S, where S is 
the three-dimensional space orthogonal to . Thus, the second fundamental form of B is 

#a/ = Kpn» - l a pv? • (3-35) 

In another basis w M , ra M for the spacetime orthogonal to B, the second fundamental form 
becomes 

= fca/jn" - 4/3^ , (3.36) 

where £ a p is the extrinsic curvature for B embedded in T (which is orthogonal to n M ), and 
k aj 3 is the extrinsic curvature for B embedded in X (which is orthogonal to u^). By using 
the boost relations (|l.l| a,b) we find that the n M and w M components of H a p p are 

k a p = jk a p + ~jvt a p , (3.37a) 
4/3 = 74/3 + ivk af3 . (|3.37| b) 

Recall that the energy and normal momentum densities are defined by e = k/n and — 
—£/ k, respectively. We therefore see that the traces of Eqs. (|3.37| a,b) yield the boost relations 
( |3.18| a,b). Also recall that the shear tensors are defined by Kr] ab = k ab — ka ab /2 and K( ab = 
_£ab _|_ 1^ j2. The trace-free parts of Eqs. ( |3.37| a,b) are then seen to yield boost relations 
(p|a,b). 

The boost invariants among e, 7] ab , and ( ab are scalars constructed from the second 
fundamental form H a ^. Note that contraction between upper and lower indices gives zero, 
since u^a^ = and n^a^ u = 0. Thus, nontrivial scalars are formed from 

Hap^H^s^ (3.38a) 
^\ P ^vH a p^H^s v ( |3.38| b) 

by contracting the free indices in various ways. For example, the invariant M 2 is obtained 
from Eq. ( 3.38 a) by contraction with while the invariant (Mi) 2 + M 2 /2 is obtained 

from Eq. ( |3.38| a) by contraction with a 011 a 136 . The invariant (M2) 2 is defined by 

(M 2 ) 4 + (M x ) 4 + M 4 /4 = H a ^H^,H l5 v H apu (5°J P p - a a pa ap ) . (3.39) 

Finally note that one may obtain (M3) 2 via contraction of Eq. ( |3.38| b) with 
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IV. CANONICAL THEORY 



In this section we consider the Hamiltonian formalism as it pertains to our bounded 



spacetime region j\A. We begin by casting the Trace-K action (|1.2|) into canonical form 
and examining the canonical variational principle. Next, we compute the variation of the 
gravitational Hamiltonian, using a lower-dimensional version of the lemma proved in Section 



II. A which was instrumental in computing the variation of the Hilbert action (2.1). Finally. 



we show that the boost relations ( 0TTS ) and ( |3.19| ) can be obtained from the canonical theory. 



A. Canonical action principle 

To write the Trace-K action (|1.2|) in terms of the canonical variables (hij, P %3 ), we first 
insert the space-time split of the spacetime curvature scalar 9ft, pl~|.|l| 

$1 = R + K^K^ - K 2 -2V fl (Ku>* + a») , (4.1) 

and the space-time split (|3.7|) of the T extrinsic curvature 0^ into the action ( |1.2|) , thereby 
finding 

S= — [ d 4 x N\fh(R + K„ V K» V -K 2 )-- ( d 3 x Ny/ZOyk - u - W) - - [* d 2 x^6 . 
2k Jm k Jr k Jb" 

(4.2) 

In deriving this expression, we have used Stokes' theorem, the result K = h^K^ = — (a lJ-u + 
rz/'n^V = £ + u ■ b along with n ■ a = ■ a, the standard identity \J—g = Ny/h, and 
Eq. (|3.21| a). The extrinsic curvature terms in Eq. (|4.2| ) can be written in terms of the 
gravitational momentum according to the relationship 

NVh[K (tt/ K" v - (K) 2 }/(2k) = P'-'h,, - 2P l W l V J - (kN/VK) [2P l 'P l3 - (P) 2 ] , (4.3) 

which reduces to an identity when the definition Eq. (|2.3| a) for P^ and the kinematical 
expression 

Kij = ~2^N {kl] ~ 2D{%V])) (44) 

are used. From the results derived in Appendix A, Eq. (|A7|a) in particular, the term involving 
the gradient of 9 can be written as 

Nu ■ ve = - a^v v )v^e = 6- v a d a e . (4.5) 

Putting these results together and performing an integration by parts on the 9 term, we find 
S = d 4 x{P ij h lj -NU- Vm % ) - J_ d*x{{6/K)^ + yft(Ne - V a j a ) } (4.6) 

where Tt and Tii are the Hamiltonian and momentum constraints, respectively. In the T 
term of S above, N = A/7, V a = V a , and e and j a are given by Eqs. ( |3.13| a,b). 
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An alternative expression for the boundary terms of S can be obtained as by using the 
kinematical equation ( J4.4p to rewrite the quantity yfa = dy/a/dt. Projecting Eq. ( |1.4| ) 
onto B, we have the result = —a^a^(a a /3 — 2D a Vp)/ (2N), whose trace yields dy/a/dt = 
y/a(—N£+a^DiVj). In the last term of this expression, a^DiVj can be simplified by splitting 
Vj into its normal and tangential parts. This yields a 13 Di((jj k V k + Nvrij) = d a V a — Nvk, 
and results in the useful expression 

dyfdjdt = y/^(KN] h + d a V a ) (4.7) 

for the time derivative of yfa. Putting these changes together, we find that the action 
equals 



S = J d 4 x{P ij h ij -NH- V*?U} - j_ d 3 xx/a{Ne - V a j a + Nj^o) . (4.8) 



In both forms (|4.6|) and (|4.8|) for the action, the independent variables are hij, P tJ , N, and 
V\ 



The variation of the action ( [4.6| ) or (4.8) can be computed explicitly, although the cal- 
culation is difficult, 12 and the result is 

5S = (terms that give the canonical equations of motion) 

rt" 1 rB" 

+ d 3 xP lJ 5h u -- d 2 X05y/^ 

Jt' K JB' 

+- / d 3 xy/^\KNj h + d a V a - (l/2)a ab & ab }S9 
k Jr 1 > 

+ J d 3 x^{-e5N + j a 5V a + (N/2)s ab 5a ab } . (4.9) 

This is not an unexpected expression, in view of Eq. ( |3.1| ) and the definitions ( |3.13|) . Notice 
that 9 need not be held fixed in the canonical variation principle, as the term which multiplies 
58 in ([4.9|) is ( |4.7| ) which vanishes as a consequence of the canonical equations of motion. We 
remark that in obtaining the result ( [4.9|) from ( [4.6|) or Q4.8|) , one must use the kinematical 
relations (|4.5|) and (|4.7|). These relations are included among the equations of motion. 



B. Variation of the Hamiltonian without boundary terms 

The "base" Hamiltonian for general relativity, unaugmented by boundary terms, is 

#base = / d 3 x(NH + V'Hi) , (4.10) 
where the Hamiltonian and momentum constraints are 

n = 4=(2p«p« - p 2 ) - ^VhR , (4.ii) 

\Jh 2k 
Hi = 2I),1>! . (4.12) 

For convenience, we write -ffbase = Hn + Hy, where e. g. we define = / s d 3 xNri. In the 
calculation of SH below, we only keep terms that give rise to boundary terms. This avoids 
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clutter in our presentation, and in any case these are the difficult terms to isolate correctly 
in the variation. 

First consider the smeared Hamiltonian constraint, denoted H^. We have 

SH N = / d 3 xNVhh ij 5R i:j , (4.13) 

where the dots denote terms that do not give rise to boundary terms. This calculation is 
nearly identical to the calculation of the variation of the Hilbert action from Section II. 
Thus, we find 

h ij SRij = D i A i , 

A* = 2h jli D k] 5h jk , 

and Eq. ( [4.13j ) becomes 

5H N = -.. + — / d 3 xVh(D i N)A i - — I d 2 xy/aNn i A i . (4.14) 
2k js 2k JdT, 



Moreover, our proof of the lemma (|2.7|) in Section II goes through unaltered for the case at 
hand (a lower dimensional setting). Therefore, we have 

mk 1 = 25k + k ij 5h i:j + di{a)8n j ) , (4.15) 

where di is the covariant derivative on <9£. Now, the first term in Eq. Q4.14 ) involves 



(D i N)A i = 2h j[l D k] [(D i N)5h jk } - 2h jli [D k] (DiN)}5h jk , (4.16) 
so that, keeping only boundary terms, we find 

SH N = ■■■ + —[ d 2 x4a{ - 2N5k - Ntf'Sha + (d;iV)<W 

2k JdT, I 

+ n i (D j N)5hij - (n i D l N)h jk 5h jk j . (4.17) 

With the substitution = + riiUj and the useful identities n % 5aij = —<7ij5n\ 5rii = 
riin^Srij, and h^ba^ = a^Sa^, one obtains 

SH N = / tfxy/al2N6k + [Nk ij + n k (D k N)a ij }5aij\ (4.18) 

2k J dT, I J 

for the contribution to 5H from the term. 

Now consider the smeared momentum constraint, denoted H^. It is straightforward to 
show that 

8Hi = -2D j [5(P> k h ki )] + P jk D t (5h jk ) , (4.19) 
from which one obtains 

5H 9 = ■■• + / d 2 x(^/Vh)^-2niV k 5P ik + [(n • V)P jk - 2 ni P ij V k }5h jk j . (4.20) 
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The first term in SH-p can be rewritten by noting that the factor y/arii/yh is metric inde- 
pendent, and can be passed inside the variation 8. With the shorthand notation 



J 



-2riiP lk /Vh 



(4.21] 



the first term in the integrand of 8 Hp becomes V k 8(y/aj k ). The remaining terms in 8Hp 
can be rewritten with the result 



5h jk = hjhfShem 



2rijn k n e 8n£ + a e j a k n 8ai m + 2a l (j n k) n m 8a, 



Cm 



(4.22) 



which is derived by using the substitution h%j = aij+riiTij and the useful identities mentioned 
previously. With these changes, we obtain 

6H 9 = --- + J d 2 x^V k 8{^ 3 k ) + v^(n • V)3 k n k n i 8n i 

+^\{n ■ V)P M aiai/Vh + (V k ai)j l aj + {V k ai)j l nM8a i3 \ . (4.23) 



Our next task is to simplify the term V k 8(y/aj k ). Using the useful identities, we find 

V k 8(V^j k ) = Vi(ai + rfn k )8{^f{o k + n e n k )) 

= (V k al)8(^j e a u ) + (n ■ V)8(^j k n k ) - ^{n ■ V) 3 k n k n i 8n i 



-^/a{V k ai)/a(8(T ij - ^(V k a i k )fn e n j 8a i;j . 



(4.24) 



The last three terms in V k 8(y/aj k ) cancel other terms in the integrand of SHp, leaving us 
with 

8Hp = --- f d 2 x\(n ■ V)(^a i k ajP H /Vh)8<j ij - 2(V e a i e )8{ v ^n J P jk (x kl /Vh) 

- 2(n ■ V)8(V^n i P ij n j /Vh)^ . (4.25) 

Here, the definition ( [4.21 ) has been used to express j k in terms of P % K 
Collecting the results from Eqs. ( |4.18| ) and ( [4.25 ), we have 



sh 



base 



+ 



/ d 2 x\ - -Ny/a5k - 2(n ■ V)8(^an i P ij n j /Vh) - 2V l ai8(y/an j P jk a ki /Vh) 

JdT, I K 



2k 



[Nk ij + n k (D k N)a i] ] - (n ■ V)alajP ke /Vh 



Saij 



(4.26) 



We can rewrite this expression in terms of coordinates x a on the surface <9£: Let <9£ cor- 
respond to an r = const surface, and define a* = dx l (x a ,r)/x a . Then 8a % a = 0. Also let 
a ij = a o,b a i^ yi = ya^ + ^ . y^i ^ and k a = k ah o l a o{. Then we have 

8H hase = ■■■+ f d 2 x\--N^8k - 2(n • V)8{y/an i P ij n j /Vh) - 2V a 8( y V°n j P jk a ka /Vh) 



-\/a 



2k 



[Nk ab + n k (D k N)a ab ] - (n ■ V)a a k a\P u / \fh 



8a, 



ab 



(4.27) 
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In terms of the boost velocity v = (n ■ V)/N and the quasilocal densities 

e = k/n , 

j h = -2niP ij nj /Vh , 
j a = -2a ai P ij nj/Vh , 



we obtain 



s ab = _ k ab + 
K \ 

t ab = 2<(x;p 4 7v / ^ 



n%N 
N 



-k 



a 



ab 



d 2 x\N 



cm 



(4.28a) 
0b) 

(He) 



(J(V?e) - v6(VZfr) - V a 5(^J a ) + (N^/2)(s ab - vt ab )8a, 



ab 



(4.29) 



for the boundary terms in the variation of the base gravitational Hamiltonian. 



C. Boost Relations for e, and j a from the Hamiltonian 

The gravitational Hamiltonian H[N, V] whose values are the quasilocal energy density e 
and quasilocal momentum density jj is 



(4.30) 



H[N, V] = d z x(NH + + J d*xy/a(Ne - V+jt 

The variation of this Hamiltonian with respect to the canonical variables is 



5H[N, V] 



J 6H[N,V] SH[N,V) 



SP ij \ + SH[N, V] 



cm 



where 



5H[N, V] 



5h 



^( P Mp ke _ p*/2)h ij + ^(2P ik P£ - PPV) 



+ 



VhN 
2k 



G 



—D { D j N + —(D k D h N)h ij 



+2(D k V (i )P j)k - D k (P ij V k 
6H[N, V] 2kN 



5PV " Vh 
and the boundary terms are 



(2^ 



Phu) + 2D (l V 3) , 



SH[N, V] 



;m 



cm 



d 2 xV^ e6N - j h 5(Nv) - j a 5V a - (N/2)(s ab - vt ab )8a ab 



(4.31) 



(4.32a) 
111b) 

(4.33) 



The terms are readily found using the results obtained in the last subsection for the variation 
of H hase . Now let us compute the change in the Hamiltonian corresponding to a quasilocal 
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boost. That is, perform a surface deformation that becomes an infinitesimal pure boost at 
the boundary <9£ (or, more precisely, becomes an infinitesimal pure boost in the orthogonal 
complement to the tangent space of each boundary point). The surface deformation is 
described by a deformation vector, which we split into a normal part (lapse function) rj and 
a tangential part (shift vector) zA The characteristics of an infinitesimal boost at <9E are 

V \ dE = , (4.34a) 

^%e = 0, (Hlb) 

and 

n%v\ d E = , ( 4 -35) 

where 9 is the velocity parameter (see Appendix C for an explanation of the relevant geome- 
try of this assignment). Under a surface deformation, the changes in the canonical variables 
are 

5h .. = htjdt = 6H 5 [ pJ ] dt , (4.36) 

5P i] ee P ij dt = - 6H }j l, ^ dt , (4.37) 
any- 
where SH/SP 1 ^ and SH/Shij are given by Eq. (|4.32|) . 

A surface deformation only affects the canonical variables. By definition, the lapse and 
shift remain unchanged, SN = = 5V l . In the surface terms of Eq. Q4.33|) , we must consider 
the variations of Nv, V a , and a a b- First, let us look at 8(Nv): 

= (n i n j n k 5h jk /2)V i \ d z 

= kn-V)n^ 6H[ \ l/] dt 
2 v ; 5PJ k 



as 

NvtfriiDjV^dtlgx . (4.38) 



In deriving this result, we have used Eqs. ( |4.32[ b) and ( f4.34| a). Now turn to the variation 



of V a . Recall that V a = V l af with cr a = a ab hijai, a ab = athijal, and with al metric 



independent. A calculation similar to the one above, which uses the conditions ( (4.34| ), 
yields 

5V a \ dJ] = Nva^n^Djiy^dtl^ . (4.39) 

Although it would seem to us not necessary, we find it convenient to choose the shift part 
of the deformation, z/, so that 5(Nv) = = 5V a . Therefore, we impose 

^'£^1^ = 0, (4.40) 

as an additional condition, along with Eq. (|4.34| ). Finally, consider the variation of a a b- It 
is not difficult to show that 



24 



(4.41) 



where d a is the covariant derivative on dH. Since v % vanishes on <9£, we find that Sa a b = 
on <9£. This, along with the results SN = 0, S(Nv) = 0, and 6V a = on <9£, implies that 
the boundary term (|4.33| ) is zero under the variation defined by the boost rj, v % . 

The results above show that the variation in the Hamiltonian is SH[N, V] = Hdt where 



H 



:i ( 5H[N, V] 8H[r], v\ _ 5H[N, V] 8H[ V , v\ 
>: I Shij JP^ 6P^ Shij 



d 3 x< 



(4.42) 



We comment on this equation in more detail below. The next step is to insert the results 
from Eq. ( |4.32j ) into the expression ( |4.42j ) for H and simplify. Although the calculation is 
essentially straightforward, it is also somewhat long and difficult. The result is 

H[N, V] = £ d 3 x^NH + Vm, 

+ / d 2 x^a\^(ND i T 1 )P ij n j + -(n l D i N)d a v a - -(rfDiNMn ■ u)k 

JdT, I y/h « « 

-(DiN)k iS i/j - -(D i N)d i (n ■ v) - -^r]{n ■ V)H 



K 



K 



n^di^crjUi) - (n ■ v)n {k k^ 
+n lk k i]e i/n j + n [k d k] (n ■ - (N, V <-> rj, v) \ , (4.43) 



+-Nn i FC i v j + JLy*pi* 



where we have defined 



N = V'D.r] - u l DiN , 

V i = ND'r) + V'Dji/ - -qD'N - uWjV 1 . 



(4.44) 
(4.45) 



Next, we impose the boundary conditions ( |4.34j ) [note, we do not need to use Eq. ( |4.40| )1, 
and obtain 



H[N,V] = J^d i x{NH + Y'H, 



+ 



d 2 x^{n i D i r j ) \2N ni P ij nj /Vh + (n • V)k/n - d a V a / '« 



(4.46) 



With the definitions ( 4.2^ ) of the quasilocal energy and momentum densities, one then writes 



H[N, V]= d 3 x\ NH + iTHi \ + / d 2 x v ^\-N6j h + Nvk + V a d a 6/K 



dY. 



(4.47) 



where # is the velocity parameter defined in Eq. (|4.35| ). 

When the constraints hold, 7i = = Tii, the energy surface density e and momentum 
surface density ji are the values of the Hamiltonian H[N,V] associated with various lapse 
functions N and shift vectors V 1 . That is, the coefficients of N and V % in the surface 
terms of H are e and ji, respectively. The expression for H above gives the change in H 
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under a surface deformation that becomes a boost at the boundary <9£. From Eq. ( f4.47| ) we 
determine the changes in the values of H under a boost, namely 



a = -9 Jh , (4.48a) 
jh = -Be , (jg b) 

j = -dj/K . (S) 



These expressions can be integrated to obtain the boost relations for a finite boost, as follows. 
The first two equations, 

ds . 4 , , 

To = ■ < 449) 



for the energy and normal momentum surface densities, have the solution 



(4.50) 



e(0) = T e(O)- 71^(0) , (4.51) 
>(0)=7>(O)-7^(O) , (4-52) 

where 7 = cosh^ and jv = sinh^. Similarly, Eq. (|4.48| b) yields 

Ja(0)=j a (0)-d a 9/K (4.53) 

for the tangential components of the momentum surface density. These results are equivalent 
to the boost relations ( |3.18| ) and ( |3.19| ). 

V. EXAMPLES OF QUASILOCAL ENERGY AND MOMENTUM 

We now examine our quasilocal energy and momentum for several spherically sym- 
metric scenarios, including static solutions of Einstein's equations, a boosted foliation of 
the Schwarzschild geometry, and isotropic cosmologies. (Dadhich and Bose have used the 
quasilocal energy and the gravitational charge defined by the Komar integral to charac- 



terize black-hole horizons for spherically symmetric spacetimes. |20|) We then examine a 
non-spherically symmetric scenario, namely cylindrical gravitational waves. These exam- 
ples, and our treatment of energy-momentum at spatial infinity in Sec. VI, require that we 
address the issue of zero-points for quasilocal energy-momentum, and we begin with a dis- 
cussion of this issue. We do not present concrete examples of quasilocal stress, however we 
draw attention to Ref. p3| in which quasilocal stress was considered by Booth and Creighton 



in their calculation of the tidal heating of Jupiter's moon Io. 



A. Subtraction term 

For any variational principle one has the freedom to add to the action terms that depend 
on the fixed boundary data. Thus, we can append a "subtraction term" — S* to the Trace-K 
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action S, which is a functional —S°[jij,hij] of the fixed boundary data 7^ and hy. The 
modified action S — S°, like S itself, yields the Einstein equations as equations of motion 
when varied subject to fixation of 7^ and hij on the spacetime boundary dM.. Certain 
modifications are brought about in redefining the action S — > S 1 — S° to include a subtraction 
term; however, before turning to the details of these modifications, let us point out that - 
as indicated in the fourth footnote of the introduction — the freedom associated with the 
subtraction term implies that our formalism alone does not completely determine a definition 
for quasilocal energy and related quantities. This ambiguity is a field-theoretic version of 
the standard one associated with any finite dimensional mechanical system described by a 
variational principle, namely the freedom both to choose the zero-point value of the energy 
and to redefine the system's momenta via canonical transformation. (Ref. jl| spells out 
the analogy between the field-theoretic freedom present in the CQF and the corresponding 
freedom in finite-dimensional mechanics in some detail. 13 ) As such, the subtraction term is 
not a background structure per se. However, we may and often do in practice introduce a 
background structure, a "reference space," as a vehicle for introducing a particular physically 
relevant subtraction term, and we therefore use the terms "reference" and "zero-point" as 
synonyms. 

The freedom to include a subtraction term in the action leads to modified expressions 
for quasilocal energy-momentum, and ones which are defined uniquely only up to reference 
contributions. In the interest of economy, let us confine our comments (for the moment) 
to the modified expression ke = k — k° for the quasilocal energy surface density, where k° 
stems from the inclusion of — S°. It should be clear that our discussion also pertains to the 
momentum densities %. The issue then is how to resolve the ambiguity in the formalism 
by selecting a suitably unique reference term k°. As we show later in Sec. VI, if the goal 
is to obtain agreement between the CQF and the accepted notions of total gravitational 
energy (either at spatial or null infinity), then there is a suitably unique choice of energy 
zero-point [|I|,|T^,|T^,^1|] . However, at the quasilocal level there is no known preferred choice, 
other than the choice k° = 0. (This simple choice has proven useful in itself, as seen 
in Sec. III.C and more recently in an approach to numerical outer boundary conditions. 
26| However, it leads to an infinite energy in large sphere limits.) At first sight it would 



seem that the the zero-point ambiguity in the CQF is no better or worse than the situation 
encountered with, say, pseudo-tensor descriptions of gravitational SEM, which are plaqued 
by ambiguity in the choice of background coordinates (or more generally in the choice of 
moving frame). However, the CQF offers some new insight into the ambiguous nature of 
gravitational SEM. Foremost, it provides a physical interpretation for the ambiguity, and one 
that is a field-theoretic generalization of the standard ambiguity present in the Hamilton- 
Jacobi description of ordinary mechanics. Moreover, as we discuss further below, in our 
formalism the selection of zero-point is usually cast as an embedding problem, affording a 
precise mathematical interpretation in terms of an associated pde problem. 

These insights aside, let us state again that we do not have an over-arching rule, ap- 
plicable for all quasilocal two-surfaces, for selecting (suitably unique) zero-points. In our 
view it is the physicist's job to select the appropriate choice of zero-point on a case-by-case 
basis, with the only guide being the rather nebulous principle that the selection should be 
tailored to the "physics" of the scenario at hand. We would like to point out that this is 
a common enough state of affairs in general relativity, a meta theory known for its wealth 
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of possible boundary conditions. Indeed, by way of analogy consider the search for solu- 
tions of the Einstein field equations. In practice, relativists certainly do not attempt to find 
the general solution, rather they attempt to find solutions given some additional physical 
input (boundary conditions, symmetries, etc.). In practice the same such additional input 
is needed to associate a meaningful QLE with a particular quasilocal two-surface. These 
considerations suggest that, rather than albatross, the zero-point ambiguity is a desirable 
feature of the CQF, as it affords pliable enough definitions of SEM to have broad application 
in general relativity. 

Turning now to the technical details, let us note that Eqs. ( |3.13j ) and (|3.15 ), along with 



the unbarred versions of Eqs. ( |3.15| ), yield the purely kinematical relationships 

k = 7A; + jv£ , (5.1a) 
I = 1 i + lv k , (|51|b) 
<t\K^ = aiK ijn j - d a 6 , (Ob) 

with similar relations stemming from Eqs. ( |3.13j c,f). The inclusion of a subtraction term, 
S — > S — S°, will modify the definition (|3.13| a) so that ke = k — k°. Here, as suggested 
in the original paper PJ, we have chosen the subtraction term S° such that the quasilocal 
energy surface density acquires a term k° which is the trace of the extrinsic curvature 
of a surface B with metric embedded in some reference space. (Note that given a 
reference space, a spacelike slice of some fixed spacetime with boundary metric equal to <j a b, 
a family of reference spaces can be generated by boosting the slice at the surface B.) This 
requires the T contribution to S° to be a linear functional of N with coefficient —^/ak°/K. 
Likewise, by choosing the S" contribution to S° to be a linear functional of M with an 
appropriate coefficient, we obtain a modified version of Eq. ( |3.13| d), namely = —1+ 1°. 
Finally, by choosing the T contribution to S° to be a linear functional of V a with an 
appropriate coefficient (and choosing the S" contribution to S° to be a linear functional of 
W a with an appropriate coefficient), we obtain a modified version of Eq. ( |3.13[ b,e), namely 
Kj a = a l a Kijfi j — {a l a Kijn^)\° . How do these modifications of e, Jj-, and j a affect the boost 



relations ( |3.18|) , ( |3.19| )? The answer depends on the relationship between the subtraction 



terms for different observers. If the subtraction terms k°, £°, and (a l a Kijn^)\° for different 



observers are chosen such that they satisfy the kinematical relationships ( |5.1| ), then the 
boost relations become 

e = - jvj h , (5.2a) 

Jh = 7Jh - ive , (CT>) 

Ja = Ja , O) 



where ke = k — k°, njy- = — I + £°, and Kj a = o l a K i j , n? — (a^KijU^l . In order to have k°, 
£°, and (a l a Kijni)\ related as in Eqs. (|5.1| ), we must choose a fiducial reference space for 
the subtraction term for some fixed observers (say, the unbarred observers), then choose the 
reference space for the subtraction term for all other observers to be boosted relative to the 
fiducial reference space. Here we do not discuss modifications of s ab and t ab arising from a 
subtraction term, but see Ref. |13j where these stresses are modified in an Ashtekar-variable 



reformulation of portions of this theory. 
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The construction of reference term k° from a reference space amounts to posing and 
solving an isometric embedding problem. One natural choice — discussed in the original 
paper Ref. - is to embed B isometrically into Euclidean three-space E 3 in order to 
obtain an extrinsic curvature tensor (k°) ab (and hence k°), a task tantamount to solving the 
following system of PDE: 

(k ) 2 - (k°) ab (k°) ab -K = (5.3a) 
d b k°-d a (k°) a b = 0, @)) 

where TZ is the B Ricci scalar and the d a denotes covariant differentiation in the B met- 
ric. These are the Gauss-Codazzi-Mainardi constraint equations for B embedded in E 3 , and 
may be derived as in Ref. |51]. Notice that {k°) ab is indeed determined solely by a ab (al- 
beit non-locally in general). With the second fundamental form (k°) ab expressed in terms 
of the embedding's coordinate chart, this is Weyl's problem, a classic problem of differen- 
tial geometry in the large for which an extensive literature exists. In a somewhat recent 
formulation of the problem, Heinz J73| has proven the existence of such an embedding if 
the B scalar curvature 1Z is everywhere positive and the metric functions a ab are of C 2 
differentiability class. Uniqueness of the embedding, up to Euclidean motions, then follows 
from the "rigidity theorem" of Cohn-Vosson. [J74] While such a Euclidean or "flat-space" 



reference proves important when one considers asymptotic limits of the quasilocal energy, 
we note that other useful prescriptions for referencing quasilocal energy can be devised; see 
for example Refs. 119,[2l|,24|. In these works, the idea has been to use a non-inertial hyper- 



plane of Minkowski spacetime as the reference space. Note that such non-inertial slices do 
not have non-vanishing extrinsic curvature tensors, and hence would give rise to reference 
contributions to the momentum densities jk- That said, in those works the main focus has 
been on the energy expression. 

In what follows, we shall assume that the quasilocal energy surface density e = (k — k°)/K 
has been set via Euclidean reference. This means that for a fixed metric geometry o ab on 
B, the total quasilocal energy 

E= — [ d 2 Xv ^(k-k°) (5.4) 

87T J B V ' 

arises as the difference between two total mean curvatures (here we set Newton's constant 
to unity so that n = 8ti). The first, the proper integral of k/(8n), is associated with 
the embedding of B in a hypersurface £ of the physical spacetime. The second, the proper 
integral of k / (877), is associated with an isometric embedding of B in an auxiliary Euclidean 
three-space E 3 (which could in turn be viewed as a slice of Minkowski spacetime). Such flat- 
space subtraction assigns that portion of the auxiliary E 3 contained within B the zero value 
of energy. Flat-space subtraction does not modify the momentum densities jk, since each of 
these vanish for a two-surface drawn in an inertial E 3 hyperplane of Minkowski spacetime. 

Finally, let us collect some overall results for referenced quasilocal energy and round 
spheres. Start with the general line-element for a spherically symmetric spacetime, 

ds 2 = -N 2 dt 2 + H 2 (dr + V r dtf + R 2 (d6 2 + sin 2 6d(j) 2 ) , (5.5) 



29 



where N, H, and the areal variable R are functions of t and r. Let £ be the interior of 
at — constant slice with two-boundary B specified by r = constant. A straightforward 
calculation of the trace k of the extrinsic curvature k a b yields 

k =-m- (5 - 6) 

where the prime denotes r partial differentiation. Now consider a round sphere with radius 
R embedded in E 3 . Such a sphere has an extrinsic curvature (fc°) a £> with trace 

fc° = -f. (5-7) 

For the scenario at hand, the referenced energy density is then 

1 /l R' 

e ~ a^\r~ 1m 

from which we find 



(5.8) 



E = R(l- R'/H) (5.9) 



for the total quasilocal energy (|5.4j). 



B. Static solutions 



Let us assume the geometry is static, with V r = and r = R. For a simple isentropic 
fluid with energy density p{R) and pressure p{R), the Hamiltonian constraint G\ = — 87rp 
implies p4 



2m 
R 



where 



H= 1-^ 



R 



-1/2 



m(R) = 4vr / dR R p(R) + M 



(5.10) 



(5.11) 



The Schwarzschild black hole solution is obtained by choosing p = p = and m = M, 
whereas a fluid star solution with p ^ must have M = for the geometry to be smooth 
at the origin. In each case, the energy is 



E = R 



1 



1 



2m{R) 
R 



1/2 



(5.12) 



with m(R) defined in Eq. ( |5.11| ). Observe that for a compact star or black hole, E — > m(oo) 
in the limit R — ► oo, which is precisely the ADM energy at infinity. p7| Section VI further 
examines the relationship between the quasilocal energy (p.4|) and the ADM energy for more 



general asymptotically flat spacetimes. We note that for the geometry at hand = 0, and 
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it follows from the results of Section III that jy- = j a = 0; that is to say, all (unreferenced) 
quasilocal momentum densities vanish. 

As discussed in Ref. |1], the Newtonian approximation for E consists in assuming m/R 
to be small, which yields 

2 

Tfl 

E^m+ — . (5.13) 
2R 

In this same approximation the first term, m(R), is just the sum of the matter energy den- 
sity plus the Newtonian gravitational potential energy associated with assembling the ball 



of fluid by bringing the individual particles together from infinity (see Ref. [54|, Box 23.1). 
The second term in Eq. Q5.13D , namely m 2 /2R, is just minus the Newtonian gravitational 
potential energy associated with building a spherical shell of radius R and mass m, by bring- 
ing the individual particles together from infinity. Thus, in the Newtonian approximation, 
the energy E has the natural interpretation as the sum of the matter energy density plus 
the potential energy associated with assembling the ball of fluid by bringing the particles 
together from the boundary of radius R. In this sense, E is the total energy of the system 
contained within the boundary, reflecting precisely the energy needed to create the particles, 
place them in the system, and arrange them in the final configuration. Any energy that may 
be expended or gained in the process of bringing the particles to the boundary of the system, 
say, from infinity, is irrelevant. 



C. Boosted foliation of Schwarzschild 



To contrast the results obtained above, let us consider quasilocal energy-momentum 
for the Schwarzschild solution, but now with respect to a (radially) boosted foliation. For 
round spheres embedded in the preferred static time slices, we have found a referenced QLE 
( f5.12j ) that is a function of the areal radius R. The boosted foliation of Schwarzschild we 
now present has an associated QLE which equals the mass parameter M for any R value. 
Such a foliation arises naturally when comparing the QLE considered here and the spinorial 
definition |60| of QLE given be Dougan and Mason. || To obtain the new foliation, start 
with the Schwarzschild line-element written in terms of the preferred static or curvature 
coordinates (T, R) , fl5"4] 



—FdT 2 + F~ l dR 2 + R 2 (d6 2 + sin 2 



(5.14) 



where F = 1 — 2M/R (and M is of course the mass parameter of the solution). In terms of 
the new time coordinate 14 



t = T + M\oi 



1-F 2 



F 2 



the line element is given by Eq. (|5.5|) with 



N = H 



V 



R, 

1(^ + 1), 
HF 2 -l). 



(5.15) 



(5.16a) 

fljrg b) 
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In Fig. 2 we plot some of the new foliation's level time slices in the right exterior region 
of the Penrose diagram for Schwarzschild. The figure depicts the t = constant slices for 




0.5 1 1.5 

FIG. 2. Penrose diagram for the exterior region of the Schwarzschild geometry. 



the time coordinate ( |5.15| ) in the exterior region of the Schwarzschild geometry. For the 
full construction and interpretation of the relevant Penrose diagram, see Ref. [^TJ. Here we 
define null coordinates (see Ref. p^| , page 832) 

U = -(R/2M - l)V2 e i?/4M e -T/4M (5 17a) 

v = (R/2M - i)V2 e H/4Af e T/4Af 5 (|5T7| b) 



where our u and v correspond respectively to w'{2M)^ 1 ^ 2 and f'(2M) -1 / 2 of Ref. ^TJ, page 
153. Moreover, in Eqs. (|5.17| a,b) the curvature coordinate T = T(t,R), as determined by 
Eq. ( |5.15| ). Then the vertical and horizontal coordinates in the figure are defined as 



P 



tan" 
tan" 



+ tan 



tan" 



(5.18a) 



Starting from the bottom of the diagram, we have drawn level-time slices corresponding to 
t = -6M, -4M, -2M, 0M, 2M, 4M, 6M, 8M, 10M, 12M. While the coordinate t provides 
a good foliation of the shown exterior region, we note that these slices cross in the past 
dynamical region. 

Now, with Eqs. ( p7L6]) and ([O]) we find K° e = = —M/R 2 for the components of 
needed to compute £ = —8iij\-, defined before in Eq. ( |3.14| a). Hence, with these components 
and with Eqs. ( |5.8| ) and ( p.!6| b) we obtain 
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8ne = %— (5.19a) 
3vrj h = ^ , (]53|b) 



again with k = 8n. It follows that both the total QLE ( |5.4|) and the total normal momentum 
J\~, the proper surface integral of both have value M for a sphere of any radius R 
embedded in one of the hypersurfaces depicted in Fig. 2. 

Now introduce a system boundary T in the Schwarzschild exterior region determined 
by fixation of the areal radius R = const, in which case the level-t spacelike slices will not 
meet T orthogonally. Observers at rest in the t slices would see a changing areal radius R, 
because with the t-foliation normal u considered as an operator u[R] 7^ 0. We can compute 
the (referenced) quasilocal energy surface density e for observers comoving with T . Although 
the energy density in Eq. ( |5.8| ) was not barred, we may nevertheless use that equation along 
with r = R and H = F 1 / 2 to find 

1 / / o i\/r\ 

(5.20) 



4nR 




the result for a radius-i? round sphere embedded in a level-T static slice. Upon proper 
integration over B, Eq. ( |5.2Up gives Eq. ( |5. 12| ) (although here we would denote the QLE 
by E). Note that because we have used —2/R both for the implicit k° in Eq. ( f5.19| a) and 
for the implicit k° in Eq. (|5.20|) , the boost relations ( [3.1S| ) are no longer valid (see the 



discussion in Sec. V.A pertaining to boosting reference terms). However, the point of this 
analysis -that our QLE (and for that matter our quasilocal momentum) is inherently observer 
dependent- should be clear. There is in no strict sense one quasilocal energy expression for 
Schwarzschild. 



D. Isotropic Cosmology 

Consider the standard isotropic cosmological model with line element 

ds 2 = -dt 2 + a 2 (t)[dr 2 + f(r)(d6 2 + sin 2 ^0 2 )] . (5.21) 

Here, f(r) = sinr, r, or sinhr depending on whether space has positive (x = +1), zero 
(x = 0)) or negative (x = —1) curvature, respectively. The quasilocal energy within a 
sphere of coordinate radius r, embedded in a homogeneous surface t = const, can be obtained 
from the general result ( |5.9| ) for spherically symmetric spacetimes. Here, we compute the 
quasilocal energy within a sphere of fixed proper area 4irR 2 (fixed areal radius R). That 
is, the history T of the boundary B is the surface a(t) f(r) = R, and the embedding 
hypersurface is chosen to be orthogonal to T. In the notation of Sec. II, we are computing 
the quasilocal energy density e as seen by the (barred) observers who are comoving with the 
system boundary a(t) f(r) = R. In the calculation below we omit the bars. 
The normal and tangent to the boundary a(t) f(r) = R are given by 

n» = (-a fa, af'/a, 0, 0) , (5.22a) 
= (af, -a fa/a, 0, 0) , (gjg b) 
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where /' = df /dr, a = da/dt, and a = (f 2 — f 2 d 2 ) x l 2 . We note that, with the definition 
for f(r), a becomes 

a = (l - R 2 (x + d 2 )/a 2 y 1/2 . (5.23) 
The mean curvature of the sphere is straightforward to compute, with the result 

k = -<r^V^n v = -\ . (5.24) 
an 

Combining this with the flat-space subtraction term k° = —2/R, we find 



E = 8irR 1 



(l - y/l - R 2 (x + d 2 )/a 2 ^j (5.25) 

for the quasilocal energy. 

The result above can be written in a more suggestive form by invoking the Einstein 
equation (% + a 2 ) /a 2 = 8np/3, where p is the proper energy density of matter including a 
contribution A/87T from the cosmological constant. The quasilocal energy then becomes 



E = 8ttR M - yjl - 8tt P R 2 /3) . (5.26) 
This result for quasilocal energy has the same form as that given in Eqs. (|5.11 , 5TT2 ) for a 



static spherically symmetric solution, provided M = and p = const in the definition ( 5.11 
for m(r). Note, however, that the matter energy density p is constant on the homogeneous 
slices t = const, but is not constant on the hypersurfaces orthogonal to the system boundary 
a(t)f(r)=R. 



E. Cylindrical Waves 



Finally, we turn to a non- spherically symmetric example, namely Einstein- Rosen cylindri- 
cal waves (see Ref. [B2] and references therein) determined by a line-element of the following 
form: 



ds 



2 _ e r-* ( 



-dt 2 + dr 2 ) + e~*rW + e*dz 2 , 



(5.27) 



where the coordinates (t, r, 8, z) range over their usual values. For the line-element (|5.27|) 
the vacuum Einstein equations become |)2 



d 2 t § - d 2 r <& - r- L d r $ = 



(5.28) 



where the function T is written in terms of the scalar field $ as 



2 Jo 



r(r,t) = - / dpp (d t $) 2 + (d p $) 



(5.29) 



We assume smoothness at the origin, so that <9 r $ — > as r — > 0. 
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Ashtekar and Varadarajan, among others, have examined the 3 + 1 vacuum Einstein 
equations under the assumption of a spacelike Killing vector field. |63|] They show via "Killing 
vector field reduction" that the system is equivalent to the 2 + 1 Einstein equations coupled 
to "fictitious" matter. Computing the total energy of the 2 + 1 system, they thereby obtain 
a 3 + 1 formula for energy per unit length along the Killing direction. The class of Einstein- 
Rosen cylindrical waves determined by the line-element Q5.27Q above -with two hypersurface 
orthogonal Killing vector fields- constitute a special case of their analysis. In this particular 
case, it is $ which plays the role of matter from the 2 + 1 viewpoint, and the expression 
they offer, 



E 



AV 



-AC 



(5.30) 



is for unit-length energy along the d/dz direction. The C in their expression is Thome's 



so-called C-energy [64], essentially just the total "flat-space energy" of the $ scalar field, 



C 



1 

16 Jo 



dr r 



(<W 2 + (d r $f 



(5.31) 



E\y can be obtained directly from the integral (|5.4| ). Consider an infinitely long cylin- 
drical two-surface B determined by fixation of r and t. B is embedded in a E hypersurface 
determined by fixation of t. Let us compute the total QLE ( |5.4|) for a section of the cylinder 
of unit proper length in the ^-direction. The mean curvature of the cylinder reads 



k 



,(*-r)/2 



(5.32) 



Now turn to the computation of the k° term. The B metric functions are easily read off 
from the B line-element, 



ds 2 = e^dz 1 + e~ *rW . 



(5.33) 



Note that over the two-surface B the field $ is a constant, and from this we easily infer that 
for the auxiliary embedding of B in E 3 the associated mean curvature of the cylinder side is 



k 



3 $/2 



(5.34) 



Hence, the energy (per unit proper length along the Killing direction d/dz) obtained from 
n is 



E 



1 



1 - e- v ' 2 



(5.35) 



Comparing the definition ([5.29 ) of T(r, t) with the C-energy Q5.31J ), we immediately see that 
the energy per unit proper length (|5.35| ) agrees with the Ashtekar- Varadarajan expression 
( f5.30| ) in the r — > oo limit. 
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VI. ENERGY-MOMENTUM AT SPATIAL INFINITY 



In this section we consider the quasilocal energy as applied to spacetimes that are asymp- 
totically flat in spacelike directions p7| , |65| -|67|], 15 and discuss its relationship with the stan- 
dard treatment of energy at spatial infinity (spi) . We begin by recalling the key observation 
from Regge and Teitelboim, |65j namely, that for asymptotically flat spacetimes the grav- 



itational Hamiltonian must have well defined functional derivatives and must preserve the 
boundary conditions on the fields. This observation leads to a two step prescription for 
building the gravitational Hamiltonian. (i) Starting with the "base" Hamiltonian (|4.10|) 
(the smeared Hamiltonian and momentum constraints), one adds boundary terms and im- 
poses boundary conditions so that the resulting Hamiltonian has well defined functional 
derivatives. For the Hamiltonian to have well defined functional derivatives the boundary 
terms in the variation of the Hamiltonian must vanish under the assumed boundary condi- 
tions, (ii) One checks that the boundary conditions are preserved under evolution by the 
Hamiltonian. If they are not, then the boundary conditions and boundary terms must be 
modified. 

In essence, the canonical quasilocal formalism is an application of step (i) above to a 
manifold with boundary. That is, by working with the action in Hamiltonian form, we 
identify the appropriate boundary terms and boundary conditions that yield a Hamiltonian 
with well defined functional derivatives. The key difference between the quasilocal analysis 
and the asymptotically flat analysis is that in the quasilocal case we do not require, as in 
step (ii), that the Hamiltonian should preserve the boundary conditions. The reason is the 
following. In the asymptotically flat case the spacelike hypersurfaces E are Cauchy surfaces. 
Thus, the data on one slice £ completely determines the future evolution of the system. For 
consistency the evolved data must obey the boundary conditions, otherwise the Hamiltonian 
on future E slices will not have well defined functional derivatives. On the other hand, in 
the quasilocal context, the surfaces E are not Cauchy surfaces. They do not carry enough 
information to determine the future evolution of the system (the spacetime interior to DM). 
Therefore, we see that in the quasilocal case step (ii) cannot, and should not, be taken. 

Note that the boundary conditions for our quasilocal Hamiltonian ( |4.3(J| ) and for the 
Regge- Teitelboim Hamiltonian for asymptotically flat spacetimes [65] both include fixation of 
the boundary two-metric, lapse function, and shift vector. In the asymptotically flat case, the 
boundary is at SPI and the boundary values are determined through the specified asymptotic 
behaviors of the spatial metric, lapse, and shift. Although the boundary conditions for the 
two Hamiltonians are not specified in the same manner, they are at least "compatible" 
with one another. Also note that with the choice of flat space subtraction, the quasilocal 
energy (obtained from our quasilocal Hamiltonian with N = 1 and V % = on the boundary) 
vanishes for flat spacetime. Likewise, the ADM energy (obtained from the Regge- Teitelboim 
Hamiltonian with N — > 1 and V 1 — > at infinity) vanishes for flat spacetime. 

To summarize, the quasilocal Hamiltonian and the Regge- Teitelboim Hamiltonian are 
both constructed from the base Hamiltonian by adding boundary terms and imposing "com- 
patible" boundary conditions, so that their functional derivatives are well defined. Further- 
more, in both cases, the energy obtained from these Hamiltonians is referenced to zero for 
flat spacetime. Given this close correspondence between the quasilocal and asymptotically 
flat analyses, it is not at all surprising that for asymptotically flat spacetimes the quasilocal 
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energy agrees with the ADM energy in the limit that the spatial boundary B is pushed 
to infinity. Although this result is not surprising, we find that a careful demonstration is 
nevertheless interesting and enlightening. 16 The remainder of this section is devoted to this 
demonstration. 



A. Asymptotic flatness 

For spacetimes that are asymptotically flat towards SPI, the E gravitational initial-data 
set (hij,K l i) obeys the following fall-off conditions: 

h ij = fij + a ij (6.1a) 
K lj = 0(r- 2 ) . (Ipb) 

Here is a flat background three-metric and an 0(r _1 ) perturbation thereof. The 
radial coordinate r = (fijX i x^) 1 ' 2 is defined in terms of coordinates x k that are Cartesian 
with respect to the metric. In this section we assume that the two-surface B is a level-r 
surface, that is to say a round sphere in the /y metric (although not quite round in hij). 
The scalar curvature of B then obeys 

U = 2r~ 2 + 3 7Zr~ 3 + 0(r~ 3 ~ e ) (6.2) 

(e small and positive). We demand that r is large enough such that 1Z > everywhere on 
B (see comments near the end of Sec. V.A pertaining to Weyl's problem). Adopting the 
usual distinction between "little-oh" and "big-oh" notation, we could write o(r~ 3 ) instead of 
0(r~ 3 ~ e ). The point being that o(r -3 ) means "decay faster than r -3 order" in this context. 
We call such a two-surface a large sphere. 



B. Momentum at SPI 

Momentum is identified with the phase space generator that moves the fields along the 
orbits of a spatial vector field £\ In moving the fields along £*, the changes in the fields at a 
given point are given by minus the Lie derivative along £\ Thus, the quasilocal momentum 
in the £ l direction is defined as the on-shell value of the Hamiltonian (|4.30|) with vanishing 
lapse N = and shift vector given by V 1 = — £ l : 

Pt = H[0,-£]= f d 2 x^C3i- (6.3) 

J B 



Here, ji = —2n k Pki/\fh is the quasilocal momentum density Q4.21p . As discussed in Sec. V.A, 



with Euclidean reference there is no reference contribution to ji, because as it stands ji would 
already vanish were S an inertial E 3 hyperplane of Minkowski spacetime. In other words, the 
reference contribution — 2(n k P ki / \fh) |° to ji vanishes for the Euclidean subtraction at hand. 
This is so because an inertial hyperplane of Minkowski spacetime has a vanishing extrinsic 
curvature tensor. Now consider the case in which space is asymptotically flat, and the 
boundary B is pushed to infinity. If the vector £ to chosen to be an asymptotic translation, 
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then we immediately find that the momentum above agrees with the ADM momentum at 
infinity. p7,p5| If the vector £ is chosen to be an asymptotic rotation, then the momentum 
P^ agrees with the angular momentum at infinity. [^5j [Although here appearing in the 
context of asymptotic flatness, note that the integral ( |6.3|) is the most general expression for 
quasilocal momentum in the £ J direction, if we take ji = —2n k P k i/\^h + 2(n k P k i/\/h)\ , now 
assuming that — 2{n k P k ij \fh) |° need not arise from Euclideain reference and hence may be 
non-zero.l 



C. Energy at spi 



Consider the quasilocal energy (|5.4j ), now with B a large sphere in the sense described 
above. 17 Let us show the equivalence between the QLE and the ADM energy. Note that B 
is a Riemannian submanifold of the Riemannian space E. Let us examine the geometry of 
this embedding in order to obtain an asymptotic expression for the integral ( |5.4|) in terms 
of the E Riemann tensor. Recall the Gauss-Codazzi-Mainardi constraints |51| which relate 
the intrinsic geometry a ab (or 1Z) and extrinsic geometry k ab (of B in S) to the E Riemann 
tensor. Among these is the following: 

TZ-k 2 + k ab k ab = a lk ^ l R ljkl . (6.4) 

As the E Riemann tensor is completely determined by the Ricci tensor Rij in three dimen- 
sions, we may replace the RHS of Eq. ( |6.4[ ) with Ricci terms using the identity 

a ik a jl Ri jkl = R- 2R hh , (6.5) 

where h denotes contraction with the B normal n k . Splitting k ab =: \ka ab + (k TF ) ab into its 
trace and trace-free pieces, we rewrite Eq. ( |6.4|) as 

K - \k 2 + (k TF ) ab (k TF ) ab = aVR ljkl . (6.6) 

From this equation we may easily obtain the relevant leading order asymptotic relationship 
between k and the E Riemann tensor. Now, the fall-off of the E metric hij determines that 

a ik a jl Ri jkl = 3 (a ik a jl R ijkl )r- 3 + o(r~ 3 ) , (6.7) 

with the superscript 3 meaning take the 0(r~ 3 ) piece of that inside the parenthesis. Now 
we make the Ansatze 18 

k = -2r~ l + 2 kr~ 2 + o(r~ 2 ) (6.8a) 
(k TF ) ab = 0(r- 2 ) (|Jb) 



Using the Ansatze ( |6.8j ) along with the expansion ( |6.2| ) in Eq. ( |6.6| ), we quickly find 

2 k = l 3 (a ik a jl R ijk i-TZ). (6.9) 

Along a similar line, we can compute the asymptotic form of k°, also appearing in Eq. ( ^.4|) . 
Indeed, striking the Riemann curvature term from the RHS of Eq. ( |6.6| ) and replacing all fc's 
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with k 0, s, we get a valid constraint for the auxiliary embedding of B in E 3 [see Eq. (|5.3| a)1. 
An asymptotic analysis of the resulting equations shows that 



2/„0 



(6.10) 



It follows that the quasilocal energy surface density e = {k — fc°)/(87r) obeys 

87T£=±VV%M) + o(r- 2 ) 



(6.11) 



in the asymptotic setting described here. 

Eq. Q6.ll ) shows that the SPI definition of energy determined by the canonical quasilocal 
energy ( |5.4j ) is equivalent to the (rather symbolic) integral 



Eoc = 77- / d 2 xy/ara tk a jl R ijkh 
16n Joo 



(6.12) 



where we use ^ as shorthand for lim^oo J B so that yfa tends to r 2 sin 9. We can rewrite this 
integral in a familiar form. Recall that ^tijki = Rijki + terms quadratic in the extrinsic tensor 
Kij [and hence 0(r~ 4 ) terms], where ^tijki is the £ projection of the spacetime Riemann 
tensor JJ^ak (the same as the Weyl tensor C^ U \ K near spi). Moreover, we may use r ~ 
^JA/(4ty) where A is the area of B. Therefore, the integral above agrees asymptotically with 
what Hayward |59 calls the Ashtekar-Hansen expression 19 



Ear = -j-\hir I d 2 xV^^ x a UK C, uXK 

87T V 167T Jb 



(6.13) 



which is quite similar to Hayward's mass Eq. (|3.32 ) . The Ashtekar-Hasen expression agrees 
asymptotically with 



£adm = ^ j> d 2 xy/an k (d 3 a kj - d k a J ^ 



(6.14) 



the familiar ADM energy. We can replace 



an 



with rx k sin 6 in this integral. 



Although apparently a result known for some time |[71|| , let us briefly argue why the 



integrals ( |6.12|) and ( |6.14|) agree in the SPI limit considered here. A more detail account of 
this issue and others concerning surface integrals at spi will appear elsewhere. [[T2|] First, 
replace the Riemann curvature term in the integral ( |6.12|) in favor of Ricci terms using the 
identity ( |6.5|) . Next, via standard techniques, obtain the identity 



Ri 



d p d ii a j) 



\d k d k a 



ij 



-dda p 

2 w i w j ^ p 



(6.15) 



where the d k s are Cartesian partial derivatives. Note that one is taking the 0(r ) piece 
of the RHS of the above equation. Finally, use this last identity along with the useful fact 



—r 



that d r a 

( |6.14| ) differ by a pure divergence on the unit sphere. 



«»j + o[r' 



) , to show that -at leading order- the integrands in ( |6.12| ) and 
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APPENDIX A: KINEMATICS OF A DOUBLE FOLIATION 

Consider a double foliation of spacetime into t = const surfaces and s = const surfaces, 
where t = const is spacelike with leaves £ and s = const is timelike with one of its leaves 
T. Let 

?V = -iW M t , (Ala) 
n„ = MV^s , (Ob) 

denote the unit-normal covectors of the t = const and s = const surfaces, respectively. Note 
that N and M serve to normalize these vectors, so u ■ u = — 1 and n ■ n = +1. The induced 
(projection) metrics for the t and s foliations are 

h/jiv = 9^v + u^Uv , (A2a) 
lnv = g^u - n^n u , 0>) 

respectively. Also define the new unit vector fields 

n M = MD^s = Mh^V v s = (M/M)h^n u , (A3a) 
u, = -NVJ = -N%V v t = {N/N)%u u , (Pb) 

where is the covariant derivative on t = const, and T>^ is the covariant derivative on 
s = const. M and N are chosen so that n ■ n = +1 and u ■ u = —1. The covectors n^, 
lie in the t = const, s = const surfaces, respectively, and are orthogonal to the two- 
surfaces formed by the intersections of the double foliation. In general and are not 
three-surface orthogonal. For example, whereas is indeed two-surface orthogonal as a 
vector field within a single s = const three-surface (as the future-pointing normal to the B 
foliation), u fl need not define a foliation of A4 into spacelike three-surfaces. The spacetime 
representation of the B two-metric may now be expressed in two ways: 

<V = 9 tu> + u ^ u v - n p. n v = 9iiv + u^u u - n^fiv . (A4) 
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Now define a time flow vector field t M along the s = const surfaces by the conditions 

= 1 , **%a = . QE^a) 
Likewise, fix a space flow vector field s M along the t = const surfaces by the conditions 



1 , s"V M t = . 
Then the shift vectors for the double foliation are given by 



The time and space flow vector fields can also be written as 



Mff + W 11 



(A6a) 



(A7a) 
(Pb) 



respectively 

Using the expressions (|Al] b) for the normal and ( |A7| a) for the time flow vector field 
we find that u^n^ = (M/N)(t^ — V^V^s = — (M/M)v, where the proper "radial" velocity 
is defined by v — V ■ n/N. In deriving this result, Eqs. (|A3|a) and (|A5|a) were used. A 
similar calculation using the expressions (|A~l|a) and ( |A7[ b) for the normal and the space 
flow vector field s» yields n^u^ = -(N/M)(s^ - W^VJ = -(N/N)(W ■ u/M). Putting 
these results together, we have (using a "•" to denote spacetime inner product) 



u ■ n 



M 
M 1 



'N [ M 



(A8) 



The normalization condition for implies 

M=(D lx sg^D v s 



.-1/2 



,-1/2 



M{n^ v n v )- l/2 = M(l + (u ■ n) 2 y l/2 
M(l + (Mv/M) 2 y l/2 . 



(A9) 



,2\- 



x l 2 . This implies, from Eq. (|A8|), that 



Solving for M, we find M = M/7 where 7 = (] 
u ■ n = —jv. A calculation similar to the one in Eq. ( |A9| ) for the normalization condition 
on gives N = N/^y. Now Eq. (A"8) shows that (W ■ u/M) = v. To summarize the results 
thus far, we have 



V-n W-u 



N M 
N _M 

7 ~ N ~ M ' 



(AlOa) 



where 7 = (1 — v 2 )~ x l 2 . 

Our next task is to express the barred unit vectors in terms of the unbarred unit vectors. 
From the definition (|A3|b) of « M , we have 
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u, = {N/N)%u v 

= (N/N)(u tl -(u-n)n IM ) 

= (1/t)(<V + 7™v) > (All) 
with a similar calculation showing that 

7v = (1/7)(tv-7W„) . (A12) 

Putting these together, we find 

% = jn^ + jvu^ , (A13a) 

= + jvn^ . ( |A13| b) 

Equivalent ly, we obtain 

»fy = 7?V ~ 7 u %t > (A14a) 

m m = 7% - 7^/i , (|A14jb) 

by inverting Eqs. ( |A13| ). 

We now derive two useful expressions, one for w M in terms of t M and the other for n M in 

terms of s M . Begin with the definition ( |A7| a) for the time flow vector field t M and write the 
shift vector as = h£V v = a£V u + JV*(V ■ n). Using the formulas ( gllj ) and (|TT|a) we 
find 

Nu>* = t^ - a%V v . (A15) 

A similar calculation starting from the definition (|A7|b) for the space flow vector field s M 
yields 

Mn" = s»- a£W u , (A16) 

where equations ( |A10| ) and (|A14| a) are used. 

The foliation of spacetime into t = const surfaces is pictured in Fig. 3. The proper time 
between t = const slices, measured orthogonal to the slices, is A" dt. In the diagram the shift 
vector V 1 points to the right, along the direction of increasing x\ so the component V i is 
positive. The proper distance between the heavy dots is 

ds 2 = -N 2 dt 2 + h^dx 1 + V i dt)(dx j + V j dt) , (A17) 

where is the metric on t = const. 

The foliation of spacetime into s = const surfaces is pictured in Fig. 4. The proper 
distance between s = const slices, measured orthogonal to the slices, is M ds. In the 
diagram the shift vector points to the past, along the direction of decreasing x l , so the 
component W l is negative. The proper distance between the heavy dots is 

ds 2 = M 2 ds 2 + %(dx* + W i ds)(dx j + W j ds) , (A18) 

where 7^ is the metric on s = const. 
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s = s = ds 

FIG. 4. 



APPENDIX B: EXTRINSIC CURVATURE "SPLITTINGS" 

Let us now derive the splitting expression ( p.5|) as well as easier splitting ( p.ll|) for K^. 
Recall that the T extrinsic curvature is defined by 0^ = — ^V a n v . With the identity 

V \u v = -K^ - u^ L a v , (Bl) 

and the definitions wy = —u ■ n and 7 = (1 — f 2 ) -1 / 2 , it is straightforward to verify that 

V^n„ = 7 V^riv - jv - jv u^a u + 7 2 u u V^f . (B2) 

From this result it follows that the projection of 0^ into B is = lk flu +'jva^a^K a /s, 

where k^ u = —a^D Q n u is the extrinsic curvature of B embedded in t = const. The identity 
( |B1~|) and the definition (|3.6| b) show that the term a^a^K a p equals the extrinsic curvature 
i^y of B embedded in a surface whose normal at B coincides with n M . Thus, we find 

<^e^ = 7 ^ + 7<, • (B3) 

Next, we note that 
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n ■ a , 



(B4) 



where we have used the Leibniz rule u u V fJ ,n u 



\{u v n v ) 



n v V,jU v 



-n u V^u u and the 



definition a u = u^X7uU u of the acceleration of u u . Finally, we compute 



-a^u 13 (7 V a n^ - jv K a p - jv u a a/3 + j 2 UpV a v) 
- 7 2 a a yV a n p + 7Va>^ a/3 + 7 2 <t?V«t; , 



(B5) 

where Eqs. ( |1 . l| a) and ( |B2| ) have been used. Application of the Leibniz rule on the first 
term on the right-hand side above yields a^u^V ' a np = —a^ripV ' a vP = a^n^K a p. Collecting 
terms, we find 



(B6) 



The results (|BJ), (0), and (g§ then give Eq. flT5p. 

The expression (|3.5| ) for 6^ still retains reference to barred quantities through the 
appearance of w^ 1 and n • a. Definitions (|1 . 1|) can be used to eliminate w M in favor of m^ 1 and 
n M . The acceleration n • a can be re-expressed with the help of Eq. (|B2|) : 



n ■ a 



-u»u v 



(i^^n v - jv - jv u^a v + y*u v V M v) 
= -7 2 w W^tv + j 3 v 3 n^n u K /lu - 7V™ • a + 7 2 w M V At w . 
The first term on the right-hand side of the last line can be re-expressed as 



u u V„n 



u^n u {K^v + u^a u ) 
jv n fi n u K l j il , — 7 n • a , 



which leads to 



n ■ a 



— jv n^n v K^ v + 7 n ■ a + 7 2 ?/ 1 V^f . 



(B7) 



(B8) 



(B9) 



The extrinsic curvature term n^n u K^ v can be rewritten using the identity ( |B 1| ) followed by 
an integration by parts; this leads to n fJj n u K fJjU — u ■ b, where b^ = n u V v n^. Collecting 
results, we find 



n ■ a 



jn-a — jvu-b + u- W 



(BIO) 



This is Eq. (|3.21 a) from the main text. 

The derivation of the splitting ( 3.1 1| ) amounts to simply projecting Kij into various 
pieces normal and tangential to B. We leave this to the reader. Let us however sketch the 
derivation of Eq. (|3.21[ b). Using the boost relations ( | . 1|) and the identity (Bl|), we find 



V a u v = jvV a n u - jK au - ju a a u + j 2 n v V a v 



(Bll) 



from which the quantity n ■ b 
[with result (ggg?)]. 



-n^n^V uU v can be expressed in terms of unbarred quantities 
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APPENDIX C: VELOCITY PARAMETER 6 AND ITS TIME RATE OF CHANGE 



Let us turn to the interpretation of 6 defined in Eq. ( J4.35 ) . We have defined 6 as such 
in order that 9 is the velocity parameter. We may verify the correctness of our assignment 
by considering boosts in flat spacetime. Let 



(Cla) 



T = x sinh t , 
X = x cosh t , 

where X, T are Minkowski coordinates and x, t are Rindler coordinates. See Fig. 5. The 
x = const 








FIG. 5. 



metric is ds 2 = —dT 2 + dX 2 = —x 2 dt 2 + dx 2 . The four-velocity U(t) along x = const 
at Rindler time t is U(t) = (l/x)(d/dt) = cosht(d/dT) + smht(d/dX). The relativistic 
gamma factor between the four-velocities at t — and t is 20 7 = —U(t) -U(0) = cosht. Now 
let t = xt denote the proper time along x = const between t = and t. The rate of change 
of r with respect to proper distance along t = const is dr /dx = t = cosh -1 7. This also can 
be expressed in terms of the gradient of the lapse function: dr/dx = J Q (n l dir]) Thus, we find 
7 = cosh ^{ftdii^dt. From the definition of the velocity parameter, we have 7 = cosh# so 
that 



e(t) 



(rjd^dt 



(C2) 



It follows that rfdiT] = 6, from which we obtain Eq. Q4.35 ). Ref. [^6| discusses these results 
in terms of degenerate foliations. 
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FOOTNOTES 



1. This is not an exhaustive list. 

2. See also the related work by Freud. 



3. For a comparison of the Trace-K action with the tetrad action of Goldberg, see Ref. [12 



4. With the Trace-K action, as with any action, we are free to subtract a function of the 
fixed boundary data. In some cases, it is convenient to fix this ambiguity by introducing 
a background space or spacetime. However, this is not a required feature of the CQF. See 
Sec. V.A for further discussion. 

5. Thus both £' and S" are, as the notation suggests, leaves of the £ foliation. We have 
investigated the more general case in which £' and £" are not leaves of S, in which case 
some portions of the boundaries of some S leaves lie in the boundary elements S" or 
However, at present we find no compelling reason to allow for such generality. 

6. By introducing a partial foliation of A4 that includes T as one leaf, we can define u M and 
n M as unit vector fields in a spacetime neighborhood of T. Note, however, that as spacetime 
vector fields, and n M are not in general hypersurface orthogonal. See Appendix A. 

7. A few words concerning terminology are in order. When we apply the variational princi- 
ple, we vary the action functional among all histories that satisfy certain specified boundary 
conditions. The histories that extremize the action under such a variation are, by definition, 
the classical histories. On the other hand, a Hamilton- J acobi variation (hj variation) of the 
action is a variation among classical histories with different boundary values. 

8. Or, indeed, useful when examining any variation involving a Ricci scalar curvature (we 
have in mind the Ricci-scalar term present in the initial-value Hamiltonian constraint). 

9. For u^, of course, e = — 1, but we shall keep track of e in order to ensure that the lemma 
also holds for n M . 

10. However, the reader should resist the temptation to identify with the future-pointing 
normal of S, as in our formalism w M need not be three-surface orthogonal. It is the case 
that the £ normal agrees with w M on the two-surface B where £ and T intersect. In fact, 
this is all that we require of S, so in effect £ represents an equivalence class of three-slices 
determined by this condition on B. 

11. In that paper, the bars were omitted. 

12. The most effective way to proceed is to use the result for the variation of the Hamiltonian, 
derived below. 

13. Nevertheless, let us present some details here. Typically, the subtraction term will arise 
as a sum of functionals (smearing integrals) over each boundary element, 



hij] = S°r[%] ~ SUhij] + S°,[h 



■tji 
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where S^„ and S%> are the same functional. The E' and E" contributions to — S° give rise to 
canonical transformation of the gravitational momentum P y ' defined in Eq. ( |2.3| a). The T 
contribution modifies IP- 7 in Eq. ( |2.3| b). Modification of EP 3 ' corresponds to the zero-point 
freedom in e, J a , and s ab . 

14. See Ref. |§ for the origins of this transformation. It may also be found by carrying out 
the null limit construction described in Ref. for the particular case of Schwarzschild. 

15. A similar analysis can be carried out for spacetimes which are asymptotically anti-de 



Sitter [58 1; and in particular, by inclusion of a lapse factor into the energy surface integral, 



one can recover the Abbott-Deser definition energy in the r — > oo limit. P JH| . [?T[ | 

16. The detailed connection between the quasilocal energy and the ADM energy has not 



been derived previously in the literature. See, however, Refs. pLITl . 



17. Be careful not to confuse the metric with the Euclidean metric of the auxiliary space 
E 3 used to obtain A; , for the metric fij lives on E. 

18. These Ansatze can be confirmed directly by computing k a b in linear approximation, a 
calculation which uses the fall-off of hij and the fact that d^r is proportional to the one-form 
normal Uk of B in E. 



19. Ashtekar and Hansen |70[ define an expression at infinity via a conformal compactifica- 
tion. This integral is the physical spacetime version of the Ashtekar-Hansen expression. 

20. The four-velocity U(t) is expressed in Minkowski coordinates so that its components 
are unchanged when U(t) is parallel transported to t — 0. 
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